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A DUAL DIFFERENTIABLE EXACT PENALTY FUNCTION IN 
FRACTIONAL PROGRAMMING 


SHRI RAM YADAV*, SHIV PRASAD AND R. N. MUKHERJEE 


Department of Applied Mathematics, Institute of Technology, Banaras Hindu 
University, Varanasi 221005 


(Received 14 September 1987) 


A dual differentiable penalty function is introduced for a fractional program- 
ming problem. Some results on local and global optimal solution of the 
problem is obtained via the penalty function method. The result thus genera- 
lize those of Han and Mangasarian! on non-linear programming. 


1. INTRODUCTION 


Han and Mangasarian’ developed a penalty function via the Wolf dual of a given 
Nonlinear programming problem (NPP) and gave some results on local and global 
optimum solution to the NPP. In the present note we generalize the results of 
Mangasarian to incorporate a fractional objective function of the form /(x)/A(x). 
Observe that in case A(x) = 1, then Han and Mangasarian’s results can be obtained 
as particular cases of our results. 


2. THE GENERAL FRACTIONAL PROGRAMMING PROBLEM 


We consider here the problem, 


_ J(x) . 
(PD rere 1G) 


Subject to g (x) <0, h(x) #0 
where f and A are functions from n-dimensional Euclidean space R" into R, and g is 
1 poe 
from R" into R”. The transformation Zz) = h(x) and z=xz, (Manas? and Schaible*”*) 


reduce the problem (P) into the following problem. 


(P’) min zat =.) 


(Z,Z9) E Se 


subject to 


«() <0.{40(2)-H} <0 





*The work supported by Senior Research Fellowship of C.S.LR. 
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on ; a 
In case f and A are convex functions, the problem (P’) is a convex nonlinear program 


ming problem. We consider the Wolfe’s dual of problem (P’) given as 


(W) max L (2, Zo» ty V) 
(Z,Zo,U,v) Ee Ratetm 


subject to 


Viasz) L (2, Zou, y= Ou>d, v>0 


where 


Z Zz 
L (2, 2054, 9) = 20 f( =) + ulg (=) + vt ( Zoh (2-)- 1). 


As defined in Han and Mangasarian! we define a Penalty function 


(Q) 6(2, 2, u,v, ¥) = L (z, 20, u,v) — & ¥ Il Veeseq) LZ (2, Zo, &, V) II? 


and consider the penalty problem : 


(0) max 6 (Z, Zo, 4, V, Y). 
(Z,Z0,U,V) ne Rnte2tm 
In case f, h and g are twice differentiable functions, the exact penalty function as 
introduced above becomes a differentiable function. We have the following result for 
Stationary points of (P), (P’), (W) and (8). In the following 


C — i. Zo) E R*), 


Theorem | (Equivalence of stationary points of (P), (P’), (W) and (6))—Let f, 
h and g be twice continuously differentiable at x. 
Then, 


(1) If (%, a) is a K-K-T point of problem (P) implies that there is a 2 € R such 
that (Z, 2, @, v) is a K-K-T point of Pee 


(2) If (2%, 2, @, 0) is is a K-K-T point of (P’) then €%, 7) is a K-K-T points of (P). 


(3) (2%, 2, a, 3) is a stationary point of (W) and yrr L (z, 20, u, v)7! exists 
=> (Z, 2, i, v) is a K-K-T point of (P’) > (%, 2, @, 2) is a stationary point 
of (6) for any y. 

(4) (Z, 2, @, Disa stationary point of (@) for y « 0 and l/y is not a eigenvalue 

of Vx L (2, Z, a, ¥) => (z, Zo, U, 2) isa K-K-T point of (P’) > (Z, 2, %, 2) 

IS a Stationary point of (W), 


In the next result we modify the notation as follows : 
(u,v) = U. Then we give the local concav 
the variables ¥ and U. 


Write (z, 2) = X and 
ity for @ (Z, Zo, u, vy) = 6 (X, U, y) in both 
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Theorem 2 (Negative semidefiniteness and definiteness of @ (.Y, U, y))—Let 

z z 
G QD = (s (2), zoh (= = ey ). Let (¥, 0) be a K-K-T point of (P’). Let f, h 
and g be continuously differentiable at Y and let vxx L (XY, U) be positive definite 


with minimum eignevalue Pp > 0. Then for . = = , (X, U) is a stationary point of 
(6) and the Hessian y? 6 (X, U, y) with respect to (X, U) is negative semi-definite. If 
in addition y > - and y G (X) has linearly independent rows, then y? 6 (Y, U, y) is 


negative definite and hence (X, U) is a strict local maximum of (6). 


Theorem 3 (Stationary point of (9) as global solutions of (P’) and (W))—Let f, 


h and g be convex and twice continuously differentiable on R”, let Y7 vy? ( AP = 
Zo 


Y > y||¥|? for all (z, 73) € Rt! and Y € R"*t and some v > 0, and let y > - . For 


every stationary point (z (y), Zo (y), u(y), v (y)) £(@), z(y), Zo (vy) are independent of 
y and z (y) = Z, Zo (vy) = 2a, where (Z, Zo) is the unique solution of (P ). 


Remarks : Incase h(x) = 1 in problem (P) we get the Theorems 1, 2 and 4 of 
Han and Mangasarian! as corolaries of our Theorems 1, 2 and 3 respectively. In the 
pattern of the proof given for the non-linear programming case in Han and Manga- 
sarian' we can prove above theorems very easily. So we omit the proof here. 
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ON A NONLINEAR INTEGRODIFFERENTIAL EQUATION 
IN BANACH SPACE 


M. A. HUuSSAIN 


Department of Mathematics and Statistics, Marathwada University 
Aurangabad 431 004 (Maharashtra) 


(Received 10 November 1986) 


The aim of the present paper is to study the existence, uniqueness and other 
Properties of the solutions of a volterra integrodifferential equation of the 
form 


u'(t) + Au(t) = f(t, u(t)) 


t $ 
+ J g(t,s,u(s), J K fs, t, u(=)] dz) ds 
to to 
u(to) = uo. 


The theory of the infinitesimal generator of Co-semigroup ina Banach space, 
a modified version of contraction mapping principle, integral inequality 
recently established by Pachpatte are the main tools in our analysis. 


1. INTRODUCTION 


In this paper we consider the volterra integrodifferential equation of the form 


u(t) + Au(t) = f(t, u (t)) 


, 


+{- u PFs, ats) f K [s, t, u (t)] dt) ds, t > f >0 
0 "9 
U (to) = uo 


wotiil) 
where —A is the infinitesimal generator of Cy-semigroup T(t), t>0 on a Banach 


space B with the norm|| . ||, KE C([R,xXR,xB.B)], gE C [Rx Ry x Bx B,B] and 
eC Thien B,B] where R, = [0, co], 


. The theory of the existence, uniqueness and other properties of the solutions of 
various special forms of (1.1) have been extensively studied by using different methods 
during the past few years! 496*7910:1415916 and the references given therein. The equation 
eae | ion of many Partial integrodifferential equations 

risé in the problems with heat-flow In material with memory, viscoelasticity 
and many other physical Phenomenon®"4+10,14. The purpose of this paper is to study 
the SSE TIES uniqueness and other properties of the solutions of (1.1). The main tools 
employed in our analysis are based on the infinitesimal generator of Co-semigroup, a 
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modified version of contraction mapping principle and the integral inequality recently 
established by Pachpatte’®, 


2. PRELIMINARIES, DEFINITIONS AND STATEMENT OF RESULTS 
Let B be a Banach space with the norm || . || and —A is the infinitesimal genera- 
tor of Co-semigroup T(t), t > 0, on a Banach space. The set of bounded linear 
operators T(t), t © R,, where R, = [0, co) isa C,-semigroup on B if 
(i) T(t +s) = T(t) T(s) = T (s) T (8), t,s>0 
(ii) 7 (0) = J (the identity operator) 
(iii) T( - ) is strongly continuous in t € R, 


(iv) ||7 (t)|| & Me’ for some M > 0 and real w and t € R, (see Martin’, 
p. 276). 


Before we state our results we give the following definitions used in our sub- 
sequent discussion. 


A continuous solution u of the integral equation 


u(t) = T(t — to) Uy + [rc — s) f(s, u (s)) ds 


ae T(r —s) fg (s,2,u (2), i Kls, E,u(E))d&)dtds  .. (2.1) 
"0 ty ‘0 


is called the mild solution of (1.1). The solution u (t) of (1.1) is said to be exponen- 


tially asymptotically stable if there exist constants M and « such that 


. |lu (t)|| S M |luoll eet), tS bo. ... (2.2) 
holds for ||u,|| sufficiently small. The solution u (t) of (1.1) is said to be uniformly 
A . . 
slowly growing if, for every « > 0 there exists a constant M, possibly depending on %, 
such that 
flu (t)|| S M |[ugl| e*'~'o, fo te. ileed) 


holds for |||] < °. 


For convenience, we list the following hypotheses used in our subsequent 
discussion. 
(H,) For allt, s € [to, 4] and x;, »: € B, for i = 1, 2, there exist nonnegative con- 
stants L;, i = 1, 2, 3, such that 


IK (t, s, x1) — K (t, $, X2)|| S Za Il — X¢l| 
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Ig (t, 5, X1, 1) — & (t, 5, X25 Y2)l] S Le [llr — ell + lr — Pall] 
If (t, 41) — F(t, X2)Il < Ls |] x1 — all. 
(H,) Forevery t' > 0 and a constant c > 0, there exist nonnegative constants 
Ly (c, t') for i = 1, 2, 3, such that 
IK (t, 8, x1) — K(t, 5, X2)ll & Li (c, t’)!| x1 — xe! 
Ilg (t, 5, x1, Ya) — & (t, 5, X25 Ya)ll 
< L, (c, t’) [Ilx1 — el] + iy — voll], 
and 
If (t, 1) — f(t, x2) S Ly (e, t')Ilxr — xl 
hold for all x1, », € B, for i = 1, 2 with ||xi|| S c, ly] S ¢ and t, s € (0, ¢’]. 
We need the following lemmas in our subsequent discussion. 


Lemma 1" (p. 28)—Let B be a Banach space. Let D be an operator which maps 
the elements of 8 into itself for which D’ is a contraction, where r is a positive integer. 
Then D has a unique fixed point in B. 


Lemma 2} (p, 1157-75)—Let x (t), a (t), b(t) and ¢ (t) be real valued nonnega- 
tive continuous functions defined on R,, for which the inequality 


x(t) << xo + j a (s) x (s) ds 
a a (s) ( b (r) x (r) dr) ds 


+ J a (s) ) b(r) ) c¢ (z) x (z) dz) dr) ds 
holds for allt © R,, where Xo iS a nonnegative constant, then 


ee on fe 2 i a(s) exp (j a(r) dr) 


xX {1 + j b(r) exp ( [b (z) + ¢ (z)Jdz) dr} ds]. 
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Theorem 1—Assume that the hypothesis (H;) is satisfied. Then for um € B the 
initial value problem (1.1) has a unique mild solution u € C ({to, «J; B) for t > t» such 


that f <t<«a, Moreover, the mapping u, — u is Lipchitz continuous from B into 
C ([to, «]; B). 


In the next theorem we give a slightly different version of Theorem |. 


Theorem 2—Let the hypothesis (7.) be satisfied then for every uy € B there is 
tmax < oc such that the initial value problem (1.1) when 4) = 0 has a unique solution 
on [0, tax). Moreover, if tmax < eo, then lim |lu (t)|] = ©. 


t max 
Remark 1: It is to be noted that Fitzgibbon* and Webb'® have studied the 
existence, uniqueness and continuation of the solutions of special form of (1.1). Here 
our method and conditions on the nonlinear functions involved in (1.1) are different 
from those used in Fitzgibbon’ and Webb’®. For our proceeding we shall assume 
that solutions of (1.1) can be continued on R, = [0, ce). 


The following three theorems deals with the boundedness, asymptotic behaviour 
and growth of the solutions of (1.1). 


Theorem 3—Assume that 


IT(¢— sl SC, (2.4) 
where c is a non-negative constant. Suppose that the functions K, g and f satisfy 

\|K (1, s, u (s))I| < Mi (5) |lu (s)|| “(2.5) 

lig (t, s, u (s), @ (s))|| S Mz (t) Ms (s) [lu ()|) + l@ ODI] ...(2.6) 

f(t, u(t) < M2 (2) lla COll (2.7) 


for all t, s € [f, 0°), u, 2 © Bwhere M, (t), M.(t) and M,(t) are real valued non- 
negative continuous functions defined on [f,, cc) such that 


PM, (s) ds < 00, | My (s)ds < 00, f My (s) ds < &. (2.8) 
1 f t 
0 0 0 


Then all solutions of (1.1) are bounded on R.. 
Theorem 4—Assume that 
IT¢ —silqC,er"™ . ...(2.9) 
for a constant 4 > 0. Suppose that the functions K, g and f satisfy 


IK (t, , u (S))\| <M, (s) ee [lu (s)II (2,10) 


lig (t, 5, u (5), @(s))|| < My (t) Ms (5) em) [liu (spl + liz (SI. 2-11) 
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f(t, u (t))I| & My (t) |lu (0)!). :..(2.12) 


for all t, s € [to, co), « > Ois a constant, M, (t), M2 (t) and M3 (t) are as defined in 
Theorem 3 and satisfy the condition (2.8). Then all the solutions of (1.1) approach to 


zero as [ > oo, 


Theorem 5—Assume that 
|T (t — s)|| < C ex" -9) a eg ED 


where C > 0 and 4 > Oare the constants. Suppose that the functions K, g and f 


satisfy 
IK (t, 5, u (Sl S My (5) e*"- Ilu (| (2.14) 
le (t, 8,4 (s), @(s))l| < Ma (s) e%-> [lu (s)|| + lla (6). (2.15) 
f(t, 5, u (t))|] < Mz (t) |Iu (t)!! ...(2.16) 


for all t, s € [fo, cc), « > Ois a constant M, (t), Mp (t) and Ms3(t) are defined in 
Theorem 3 and satisfy the condition (2.8). Then all the solutions of (1.1) are slowly 
growing. 


Remark 2; Itis important to note that in Pachpatte'! has studied the stability, 
boundedness and asymptotic behaviour of the solution of (1.1) when A = 4 (t), A (t) 
is a One parameter family of closed linear Operators. Here our results in Theorem 5 
are obtained by using different approach from those in Pachpatte!!, 


3. PROOFS OF THEOREMS 1 AND 2 
In the space C = C ([fo, a], B), we define the norm 
llullc = max |lu (r)I). 
eee lu (t)|| ret fs SY a 


It is easily seen that C with the norm defined in (3.1) isa Banach space, we define a 
mapping 


FC => 'C 
by 


(Fu) (t) = TG — to) up + j T(t — s) f(s, u (s)) ds 


+{TG - JI (s,-, u(t), § Ks, E, u (E)) dé dz ds 


bo St<a, (3.2) 
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Clearly, the solution of the equation (1.1) isa fixed point of the operator equation 
Fu=u. Letu,v € C, from (3.2), (3.1) and hypothesis (H,), we get 


(Fu) (t) — (Fr) (t)|) 


< { MLs lu (s) — v(s)| ds 


4 j M § Lalu () — » (x) + fh, lu () — v || dé] de ds 


(t — %)? 
< ML; |lu — vile (t — to) + ML, |lu — i me a 


— 3 
+ ML, Li lu — vile KS 


=a Ut —1)[ Le + 1, £59 41.1, | — Yc 
A ey 


where M is a bound of ||T (t — s)|| on [t, @]. Using (3.2), (M1) (3.3) and iteration it 
follows easily that 


(Ft) (.) — (Fv) (II 
ee tM | Ls fo Ts an Bee oe A |} lu — vilo. 
...(3.4) 
From (3.4) and (3.1), we obtain 


Fru — Fvilc < Alu — vile 
where 








1 n 
A = nl | Ma | Ls +L: at hls ale 


For n large enough A < 1. Thus, there exists a positive integer n such that F" is a 
contraction in C. It follows from Lemma | that the operator F has a unique fixed 
point in C. This fixed point is the desired mild solution of (1.1). 


Let vy be another mild solution of (1.1) with v(t) = % on [fo, ]. From (3.2) 
and hypothesis (H,), we have 


lw (t) — v (I 
<M |lup — voll + § ML llu (s) — v(s)I] ds 


0 


+f f ML [lw () — ¥GI 


‘oo (equation continued on p. 522) 
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+f Liu —y Ql délds ds 


0 


where L = max {Lj, L2, Z3}. An application of Lemma 2 yields. 


lu (t) — v Cll 
< M | — vol [1 + J ML exp (j ML dz) 
x {1 + fexp (f [1 + L]dé) ds} ds]. A355) 


From (3.5) and (3.1) it follows, 


lu — vile 


‘ s 
< M |luo — voll {1 + J ML exp (J ML dz) 
t t 


0 0 


x {14+ fexp( J [1 + L] d8) ds} ds] Wen 


0 0 
which yields both the uniqueness of u and Lipchitz continuity of the mapping uy — u. 


This completes the proof of Theorem !. 


First we shall show that for every tf) > 0 and uy) € B, equation (1.1) under the 
assumptions of theorem has a unique solution wu on an interval [f9, t,] whose length is 
bounded below by 


8 (to, ||voll) 


I|tto|| 


[k (to) {La (k (to), to + 1) +°L (ke (to), ty + 1) 





= min | ie 


+ TRUK. tal) cx (Ea es 
+ (N, (to) + Ne (t))] | 3.7) 


where Lj (c, ¢), i = 1, 2, 3 are local Lipchitz constants defined as in (H,), and 
M (to) = max {|T(¢)|:0< tS 4+ 1} 
k (%) = 2 ||teq|| M (fo) 


Ny (to) = max {|| f(t, OD] :0<¢t< to + 1} 
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N, (fo) = max {llg (t, S, 0, Ik Is, 7,0 | dt)|| | eu < s < t < to + 1}. 
‘ 
0 


...(3.8) 
Indeed, Let t) = f + 8 (fo, ||oll), where 5 (fo, ||uoll) is defined by (3.7). It is a easy to 


see that mapping F defined by (3.2) maps the ball of radius k (to) centred at 0 of 
C ([to, ¢,] ; B) into itself. This is clear from the following estimate. From (3.2), (3.7) 
and (3.8), we have 


\\(Fu) (#)}| 


< M (th) |lvoll + i IT (t — s)il f(s, u(s) — f(s, Olds 
4 j IT (¢ — sh F lig (ss =, w (=), SK Es, & wD] ay 
0 0 0 


— 2 (5,7, 0, j k [x, E, 0] d&)|| de ds 


0 


‘i j IT (t — s\il f(s, Ol ds 


+ FIT (=) Fllg(, +0, fk Es, & 0} d8)| ds ds 


0 0 0 


< M (to) lluoll + M (ta) Ly (K (ta), to + 1) k (to) (t — ty) 


+ M (te) Ly (k (to), to + 1) k (to) oA Seetilc >i to) 


+ M (ty) La (h (tos to + 1) La (toate + 1) (to) SS 


< M (to) loll + M (te) (t — to) Ue (ta) (La (K (to), to + 1) 
4 Ly hte), ty +1) + Li(h (to), to + 1) Le (K (to), to + 1D} 
+ (NM, (to) + Ns (t0))] 

<M (to) |ltol| + M (to) lel 


= 2M (ty) ljuell = k (to) 


where the last inequality follows from the definition of t,. The operator ¥ satisfies a 


uniform Lipchitz condition with the constant 
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L = M(t) [Lalk (te), to + 1) +. Lak tos fe + 1) 


+L, (eto) to + 1) Le (k (tos fo + DI 


and thus by Theorem 1 it has a unique fixed point u in the ball which is the desired 
solution of (1.1) on [fo, t,]. | From above it is clear that if u is a mild solution of (1.1) 
when fy = 0 on the interval [0, +] it can be extended to [0,7 +9] with 5 > 0 by defining 
on [t, t + 3], u(t) = w(t), where w (ft) is a solution of the integral equation 


w(t) = T(t — t) u(t) + j T(t — s) f(s, w (s)) ds 


t s & 
+ f§T(t—s) J g(s,& w(E), SKE, 0, w (n) | dn) dé ds 
ee fla eee 3.9) 


A . 
Also, 6 depends only on ||u (t)||, & (<), Ni (t) and N, (+). Let [0, tmax) be the maximal 
interval of existence of the mild solution u of (1.1) when tp = 0. If thax < o¢ then 
lim lu (t)|| = co, since otherwise there is an increasing sequence f, converges to 


fe each that |/u (f,)|| < C for alln. From the above proof it would imply that for 
each f,, near enough to tmax, definition of u can be extended from [0, f,] to [0, t, + 8] 
where 5 > 0 is independent of t, and hencewu can be extended beyond tmax which 
contradicts the definition of tmax. 


For the uniqueness of the local mild solution u of (1.1) when fp = 0 we suppose 
that if visa another mild solution of (1.1) then on every closed interval of LOS aaa p 
both u and vy exist they coincide by uniqueness argument given inthe proof of 
Theorem 1. Therefore, both u and v have the same fmex and on [0, tmax), U = Y. 
This completes the proof of Theorem 2. 


4. Proors oF THEOREMS 3-5 
From (2.1), (2.4) — (2.7) we have 


lw (OI < IT — t,) nll +) IT (t — s) IS (s, « (s))Il ds 


+ JIT C= SMF lig +, uC), fk Es, & u(® | dB)ll de ds 


0 


< Cll + 5 Cy My (s) Ite (5) ds 


t 
+IC 


0 


J Mz (s) My (x) [Ie (=) + JM, @ lu @ dé] de ds. 


-».(4.1) 
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Using Lemma 2, we get 


le (t)I] < Cy fuoll | 1 + j C, M, (5) exp (J eM, (=) ds 
0 


x {1 + j Ms (=) exp ( J [Ms (6) + Mi (6)] dé) dz} ds. . (4.2) 


Thus in view of condition (2.8), boundedness of u(t) follows. This completes the 
proof of Theorem 3. 


From (2.1), (2.9) — (2.12), we obtain 


A t A 
lle (t)] S Cy e-*Fo) lull + J C, ee) Mz (5) lu (s)]| ds 
t 
0 


i A s A 
+ § C, e~*(t-9) { M, (s) Ms (t) e-* 8-9 
t t 


0 0 


(lu (x)\| +} M, (@) e-®-® |lu EI] dE] de ds. (4.3) 


Multiplying both sides of (4.3) by e*‘, we have 


lu (t)I| e# < C, eo |luoll + J C, &* M, (s) |lu (s)|| ds 


0 


t A Ss. A A 
+ § C, e%* { M, (s) Mg (+) |lu (*)l] em . e* dz ds 
‘9 "9 
t A s yey es, Ws A A 
+ fC, ew* § Mz(s) My (t) e7** e*" { My (&) e~" e*® |u(8)I dédx ds 
i t ‘o 


0 0 


A t A 

= C, e%%0 |luol] + J Cy M2 (s) llu (s)I] e*? as 
t 
0 


+ | C, My(s) | Me) [lu er teas 


‘9 0 


ECM) EM) F MO lu Ol ert ab de ds. 44) 


t 
0 


Applying Lemma 2 to (4.4), we get 
A ‘ 3 
\\u (t)I| et < C, e*"o |lvoll (1 + J c M,(s) exp ( ie M, (t) a) 


0 
(equation continued on p. 526) 
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(1 + f My (t) exp (J (Ma (E) + M, (€)] dé) de} ds], 
(4.5) 
In view of (2.8), we have 


lle (t)|| e <M, (4.6) 


where M, > Oisaconstant. Thus as t — oe, the solutions of (1.1), approach to 
zero. This completes the proof of Theorem 4. 


From (2.1), (2.13) — (2.16) we obtain, 


lu (ON) < Cie" full +f Cie My (s) lla (s)le (59) ds 


0 


‘ A s 
+ J Cy eX) | My (s) Ms (x) ex" 
t oy 


0 


X [llu(s)|] ++ J My (E) e®(-® IIe (EI dE] de ds. ...(4.7) 
0 
Multiplying both the sides of (4.7) by e-8, we have, 


\luu (t)]| em gC, e-%y I|tuol| + i C, M, (s)\) e-™ ds 


tt) 


h f C, Mz (s) f Ms (+) lu (x)|] e-** de ds 


t s r A 
+ J CM, (s) J Ms (t) J M, (€) |lu (&)\| e-®5 dé de ds. -+-(4.8) 
0 0 "9 
An application of Lemma 2, yields 


le (Me! Cr em" foul 1 + F Coat, @) exp (FAL (=) d) 


0 


(+ J Ma (=) exp (FIM, © + a, @1d8) de) do 
; , (4.9) 


In view of the condition (2.8), we obtain 


A 
Ilu (¢)|| ems < Ms. 
where M, > O is a constant, 
Ms small enough, and the pro 


...(4.10) 


The above estimate yields the desired result if we choose 
of of Theorem 5 jis complete. 
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5. EXAMPLE 


In this section we give an example to illustrate the applications of our results. 
Consider the partial integrodifferential equation of the form: 


us (t, X) = uxx (t, x) + J a(t, 5, Uxx (Ss, x), { 0; (5, t, Uxx (t, x)) dt) ds 
0 


+ a, (ft, uxx (t, x)), Oe 16> 0 nade) 
with the given initial and boundary conditions 
u (0, t) = u (1, 0) 


8 Ta ae | me Bes 
u (x, 0) = up (x) 


wheres: R, xX Ry X RX R>R, 9,:R, XxX Ry X R>R and o,: Ri X RR 
are continuous, continuously differentiable with respect to the first argument, uniformly 
Lipchitz continuous in s,+ and ft € [0, co) respectively and satisfy the following 
conditions; 


(a;) | 0, (t, s, Z;) — 91 (t, 5, Zr) | & Zi | 21 — 22 | 
(a,) 1-9 (fo18. tan 25) — ¢ (E78, 21, %). |. Lab 2. — 21) + bea 
(as) | 92 (t, 21) — 92 (t, 22) | < Ls | 2, — 22 | 

for 21, Z., Z1, Z. © Rand constants LZ, > 0, i = 1, 2,. 


We first reduce the problem (5.1)—(5.2) to the form (1.1) by making suitable 
choices of A, f, Kandg and illustrate the hypotheses of main results established in 
section 3. 


For, Let B = L? (0, 1, R); As = 2’; 
DiAjH] Ze Bz eB; z (0) = z (1) = 9} 


and K: R,x Ry X BoB; g: Ry X R, xX Bx B>B, f: R. xX B> Bare 
defined as follows 


f(t, 2) (x) = 42 (t, 2” (x), (62) E Ry xX B 
K (t, 5, z) (x) me OF (t, 5, * (x)) 
g(t, s, 2, &) (x) = 9 (t, 5,2" (*), 8), O< x <I. ..(5.3) 


Equation (1.1) then becomes the abstract formulation of (5.1)-(5.2). Now from 
(a,)-(as) and (5.3) all the assumptions of Theorem | are satisfied and there exists a 


unique solution of (5.1)-(5 2). 


Suppose that the functions involved in (5.3) satisfy the following conditions; 


(b;) | 9, (t, 8, 21) | < My (s) | 21 | 
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(b2) | o (t, S,Z1, 23) |< M, (t) Ms(s)[ | 2; | = E | Z, | ] 
(bs) | o,(t, 21) | < M, (t) | z, | 


and also 
(4) | TH -5s)|<CG, 


where M, (t), i = 1, 2,3 are non-negative real-valued continuous functions and satis- 
fying condition (2.8) and C; > 0, is a constant. Then by Theorem 3 all the solutions 
of (5.1)-(5.2) are bounded. 


Further if the functions in (5.3) satisfy the conditions; 
A 
(ci) | u(t, 5,2,)| << My, (s)'exe0 9) fa 


(ca) | 9 (t, 8, 215 2,) | < Ma (t) Ms (s) e*-) [| 21 | + | za | 
(cs) | @ (1, 21) | < Mz (t) | 2 | 

and also 
(c) [T(t —s)| < Cyemi-n 


where function M, (t), i = 1, 2, 3 are same as defined above. Then by Theorem 4 all 
the solutions (+.1)-(5.2), approach to zero as t > oo. 
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P. G. Department of Mathematics, Ranchi University, Ranchi 834008 (Bihar) 
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In this paper it is shown that the theory of linearization of a quadratic trans- 
formation in a genetic algebra Ay, initiated by Holgate® and developed by 
Abraham'”, can be extended to include some useful special cases. A general 
method of constructing basis monomials of the induced linear space B, is 
described. . Finally proceeding parallel to the original theory, the maximum 
possible dim B, as the number of solutions of a diophantine inequality and 
equality, are obtained. 


1. INTRODUCTION 


The works of Etherington‘, Schafer? and Gonshor® led to the introduction of a 
genetic algebra (GA) in population genetics. A genetic algebra is the algebra A, which 
admits a (canonical) basis {co, ¢,, C2,..., C,} obeying 


CrCy = & Nyx Ce, i,j = 0, 1, 2, 0.5.0 
k=0 


where jj, are field elements with the properties 
G (i) Aovoo = | 
G (ii) Agje = 0, for k <j 


G (iil) Ni jk a 0, for k <S max (i, j), if i, J > 0. 


. 


In a genetic algebra the transformation @:X — x® is basic representing the transition 


from a distribution x of genotypes to the filial distribution x? 


under random mating. 
Hol gate® 


made an important contribution by showing that ¢ can be linearised i.e. the 
underlying vector space A of a genetic algebra with a canonical basis and with an 
idempotent can be embedded into a higher dimensional space B (of minimum dimen- 
sion) such that ¢ has the image é (a linear transforma 
arization of a quadratic transformation in a_ genetic 
satisfactory form by Abraham!-*, Abraham! has 


(a) proved th 


tion) in B. The theory of line- 
algebra was brought to a more 


at dim B is independent of a basis of the genetic algebra 4, 


(b) described a method to generate the monomials required to form a basis of B, 
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(c) obtained dim B, in terms of dim A exactly, recursively and asymptotically. 


He has obtained, as remarked by him, the maximum value of dim B for a genetic 
algebra of order (nm + 1) by assuming that none of the multiplication constants of 
G.A.is zero. Dim Bis too great even for GA’s of smaller dimensions, e.g. for 
a= 7 + 1 = 8, dim B = 27337, vide Abraham’. However, in many concrete 
genetic situations, the multiplication constants vanish in a definite and simple pattern. 
In these cases the upper limit on dim B is greatly reduced which cannot be deduced 
from Abraham’s results. There are genetic algebras A, satisfying 


Ch = Co sala) 


and 
Aijk = 0 if 0.5, f7 a: | NY Ras RECS saUlel 


(or more generally Ay, = 0,1f0 <i,j< kandk=4t,%, ...,t) .. de) 


(1; ) expresses the fact that c; cy has no component with respect to cx if {i, 7} # {0, k}, 
k => th, Levers Pe 


For infinite random mating populations of diploids with Mendelian segregation 
such algebras are 


(i) gametic and zygotic algebras for a single locus, 
(ii) gametic algebras with uncoupled or strictly coupled loci, 


(iii) gametic algebras with coupled loci and the linkage distribution p(U) 
satisfying the condition 2 p (U) < 1. 


For the description of (i), (ii) and (iii) we refer to W6rz-Busekros*. It is interesting to 
note that if (1s) [or (1; )] is supplemented by the condition Agxx = 1/2, kK = 1, 2,...,1r 
(or k = th, to,..., tr), which is usually the case when we assume Mendelian segregation, 


then A, has an r-parametric family of idempotents [vide Theorems 4.6 and 4.8 in 
Worz-Busekros®]. In a genetic algebra idempotents represent equilibria of a population. 


The purpose of this paper is to extend Abraham’s theory to obtain a (minimal) 
linearization even when some multiplication constants of GA are zero. In Section 2.2 
and section 3 we use our method to obtain a basis and the dimension of the induced 
linear space B with respect to a quadratic transformation in 4, satisfying (1a) and (1s). 
We have stated a theorem describing two methods of generating basis monomials neces- 
sary to linearize a quadratic transformation in A,. This theorem is a generalization 
of Theorem 2 in Abraham’. However, we have followed a different line of reasoning 
for the proof. Finally we have obtained the upper limit on dim B for genetic algebras 
A, of dimension n + | as the number of solutions of a diophantine inequality and 


equality. 
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2. LINEARIZATION OF A QUADRATIC TRANSFORMATION 
2.1. Linearizing a Quadratic Transformation in a Baric Algebra 
It is known that the class of baric algebras is wider than that of genetic algebras 


[vide. Wé6rz-Busekros‘]. Hence we describe first the linearization of map ¢ : x > x? 
for a baric algebra 4,,. 


Let x = (1, x1, X%,...,%,) be an element of unit weight in A,. ¢ maps 
(1, 1, X2,..., X,) into (1 ¢, x, , X, ¢,..., Xad) where 1¢ = 1 and x,¢ are linear combi- 
nations of monomials in x,, X:,..., X, of degree < 2. Let +: be the linear mapping 


: 1 2 n 
extended to the linear space spanned by the monomials x,’ x,” ... x,, such that 


xo = Xd 
and 


4 i 


(4d x20 52) dim Ged) ce eee 


If there exists a least integer s > 1 for which the set {xid, Xk $%,..., Xk 5} is linearly 
dependent, then it spans a subspace Lx,. Let uscall it a cyclic subspace (of dim 


Ss — 1) generated by x, with respect to ¢. If the finite dimensional cyclic spaces Lx, 
exist fork = 1, 2, ... .n; ¢ is said to be linearisable, Lx + Ly +...4+ [x = B, is 


called the induced linear space. The restriction of ¢ to B, is the representation of dg. 


The set of all basis monomials of Lx,, k = 1, 2, ....m; form a basis of B,. 


2.2, Linearization of a Quadratic T; ransformation in a Genetic Algebra 


Let 4, be a genetic algebra over the real or complex field with a canonical basis 
{co, Ci, Cs, ..., C,} satisfying the conditions (1a) and (14). 


Ifx =¢9 + x, ot... +Xn C, is any element (of unit weight) of A,, then x? =c, + ne 


Qt... $x) o 


where 


g ™ QAn XK Ag = Aga: Km FLY ene 3 


and 


k-1 
ko = xR AR + ae X1 Xj Aijxe Where xy = 


= 1, =F 1) oy 
Ak => Aook 


(equation continued on p. 533) 
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k-1 k-1 
= 2Xx Ax ~ 22 Aiok xi + >> At jk XI xj. 
i—1 jm) 
Let us introduce 


l, = a linear combination of x, X., ..., Xr 


1’ = a linear combination of all monomials x; UC A Pe Aaa 
k ‘ ‘ s ‘ 5 
| = a linear combination of all monomials in Xi, X»,..., xr of degree k. 


We assume that the coefficients in /* are immaterial. Thus they obey the ordinary 


law of indices e.g. /™ 1” = 1""",(J")" = 1" and alsoa/, = 1,, a being any scalar. 
In order to get maximum possible dim B, we further assume that the coefficients of 
every lk are non-zero. 


Then 
} 
xia = 2 Ak Xk, when k = I EO a 
\ <-42) 
=f? +h,whenr<k<n ; 
2 [2 
Now ¢ : x — x? maps (1, Xj, %2, ..., ¥ ,) into (1, xt a eee ey, where x, } are 
given by (2). 
Thus 
ld = ' 
Xero = 2An xn, kK = 1, 2,..., 17 
| 
X41? = is + Irs, | 
ee i} 
X42 = apes + Thos f ( 
| 
| 
| 
x,¢2 Fatt } 

Let ¢ bea linear mapping extended to the linear space spanned by the monomials 
' — =x 
= x? Put x! " (where i’s are non-negative integers and c 1), such that xx ry Kh 

and 


6 1 


Got lt DB = (1 DS A) C2 
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Obviously the spaces Lx, are one-dimensional cyclic spaces for k = 1, 2, 


The following theorem gives the linearization of ¢. 
Theorem 1\—For k = r + 1, r + 2,..., n, Xa $** generates the basis monomials 


of the cyclic spaces Lx, which form the ascending chain 


Lx Le C; see a Lx . 
Proor : By (3) fork = 1, 2,...,r 

Xe Go = 2 Xx Xx 

> XG = 2 Ag XE 
>id = I. Pee fe 

Also by (3), for k > r, 
baa =P 4 met ...(4i) 
haod a gee + bya == Xreo d ...(ili) 


(2 — r + 1), 


Also 


Hence 
xiig¢ = igi bee (say) 
Xr42 g? = (hr ¢) +his o 
me (FF + Iea:)? + Ira + lees, by (ii) and (iii) 


= (I; i Iai)? + leas 


7 iat t+ bry, = Lrys (say) 
(equation continued on p. 535) 


au ¥Me 
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Xres ¢! — Teeacd + Iras | Liss (say) 


tote EO. +d, = L, (say). 


_ 


Also it is easy to verify that ¢ “saturates” xr4, in the sense that Xr+, pstl = 1T,..¢ 
= L,,,fors=1,2,..,n —r. Lr4, contains all the monomials which form a basis 
of the cyclic space Ls,,,° Hence we say that x,,, $f generates a basis of the cyclic 
space Lx s=1,2,..,n—r. Clearly Z,4,,1 contains all the monomials of Zr,, 


i.e. rx, is a subspace of Dy for s = 1,...,2, n— r—1, which is expressed as 1 


Lx, Cc Lx... Cc eee E> Lx . 

Consequences of the theorem : 

(a) The induced linear space corresponding to (r + s + 1) —dim genetic algebra 4,4, 
is By, = Lx 5 = 1,2,...,n—r. 


(b) If ss = x, + X2 +...+ Xr, the basis monomials of B-4,, Brio,---, B, (except 1) are 
generated successively by 


Pray (%, X2,--5 Xr) = he + Sro1 == Pry (say) 


2 
Prae (X15 Xe9-+-9 Xrto) = Pray + Sete = Pree (say) 


...(4) 


Pe (Xp Xses Fe) Poy + ha Po (say) 


ae 
ny 


or equivalently the monomials of Bris, Bry2y---, B, (including unity) are generated 
successively by 


Praia (1, *1, heey = US Sr)? + Xre1 = Pro (say) } 


Pra qi, x), Moyers Xr+2) oa Pr. To ars Prso (say) 


ya 
~~ 
un 
~~" 


PME i, Secora Po tik. (say) 


when r = 1, this mode of generating monomials reduces to one given by Abraham’. 


3. Tre DIMENSION OF INDUCED LINEAR SpAcE B, 


In Propositions 2 (a) and 2(b) below we show that dim B, can be expressed as the 
number of solutions of a linear diophantine inequality which can be converted into 


an equality. 
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a a a 
n 


: 1 2 
Propositions 2a—Dim B, is the number of monomials x, x,”... x,” such that 


@,, 4, ...,4, are integral non-negative solutions of the inequality 


ay, + as + pee ctanar + 244.1 + 2? Aras + . Tt Feds a, = j sash Ynsr 
6) 


2(b) Dim B, is the number of integral non-negative solutions of the linear 
diophantine equation: 
Ay + @, +...+ Org, + 2Arge +H oe + 2°77 Gag = 29-7, Yon > (r + 1). 
tO 
PROOF OF PROPOSITION 2 (a): We define the weight w (x;) of x;, i = 1, 2,..., mas 
w(x) = 1, fori = 1,2,..,r 
= 2‘, fori = (r+ 1),(r + Ze 
and 


¢, 4, a n 
W(x, xX, ... xX") = ZS aw(x) 
, i= 


when n = r, (6) is a, + a, +...+ a <1 which gives (r + 1) solutions equal to the 
number of monomials 1, x, X, ..-, Xr which span B,. Similarly we can show that the 
Proposition holds form = (r + 1). Next let us assume that Proposition 2 (a) holds 
for Am.1,m > r. By (5) the basis monomials of the induced linear space B,, is gene- 


rated by P,, = P*_. + xm. By inductive hypothesis P,; is spanned by all the 
monomials in the set 


a 
1 


Mmn-, a {Pa = x, sa 2 ae ee :w (Pa) < 2m-r=1}. 


[since w (pa) = 2M + ary IMG, q 207-1), 


P,, is spanned by the products Pa Po (Pa, Po © M, -1) and x,,. 
Also 


W (Pa Po) = w (Pa) + w (po) < Qm-r-1 + Qm-r-1 _ Qm-r 
Ww (Xm) = 2”—r. 


Hence P,, is spanned by all 


monomial p such that W(p) < 2™-r 
Proposition 2 (a) holds for Ay. 


By induction 
Hence the proposition, 


PROOF OF ProposiTIONn 2 (b); 


| Let D (r, n, 2"-7) denote the number of solutions 
of the equation 
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Q, + a2 +...+ ay + a4, + 2ay,, +...4 204 AE EY ala 


It is known’ that D (r, n, 2"-r) is the coefficient of t2""" in the expansion of 


én () =(1- ty — Ae) a ey 
in powers of ¢. We have 


P(r+1)s(n+1) (t) =. (1 as t) : Prsn (t) 


= (1+t+#?+...)[D(r, n, 0) + D(r,n,1) t+...4+D(r, n, 2"-") erred, 


Equating the coefficients of ¢?” " from both sides, we get 


n—r 


D(r + 1,0 +1, 2”-") - D(r, n, k). 
-0 


The right-hand side is the number of all solutions of the inequality (6) whereas the 
left-hand side is the number of all solutions of the equation (7). Hence the result. 


Remark : When r = 1, A, is a genetic algebra with restriction c2 = cy. Hence 
Propositions 2 (a) and 2 (b) obtained here are extensions of the corresponding results 
of Abraham [vide Propositions | and 2 of Art. 3 in Abraham]. 


4. FurRTHER EXTENSIONS 


(i) There are certain genetic situations in which (1; ) holds for k = t,, ty, . ,tr. 
For instance, in Sec. 1, (iii) is such a situation. In such cases the method in section 
2 can be used to show that maxium value of dim B for a genetic algebra of dim (n +1) 
is the number of solutions of the diophantine inequality. 


1-1 ir 
a, + 2a, +... + Ie (a,,-1 + iyo ee (6b -1. “+ a, ) 


iw at, ey eT. 


(ii) The method in section 2 can be used to obtain the induced linear space B, 
with respect to a (minimal) linearization of a transformation v, defined by - (x)=Go+ a 
x + ... + a, x’, where x’, where x’ is anrth principal power ina genetic algebra. 
We have established in Singh and Singh’ that B, is independent of r and coincides with 
the induced linear space B with respect toa quadratic transformation in a genetic 
algebra with non-zero multiplication constants. 
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TOPOLOGICAL GROUPS 
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(Received 17 September 1986; after revision 2 February 1988) 


Let G be a multiplicative group. Equation 

f(xy) + f(xy) = 2f @)S(y) for x,y EG +(A) 
is called D’Alembert’s Funetional Equation on G. It was completely solved 
by Cauchy for a real-valued continuous function fon the additive group R 
of real numbers with the usual topology (Aczel, 1966, pp. 117-20). In this 
paper, we solve (A) when fis a continuous complex-valued function on the 
direct product of a family of subgroups of R”, with the product topology. 


1. INTRODUCTION 


Kannappan‘ considered complex-valued functions on arbitrary G satisfying 


Sf (xy) + f xy) = 2f (x) f Y) for all x, y in G -(A) 
and an additional condition 
f (xyz) = f (xzy) for all x, y, z in G. ...(B) 


In section 3, we prove some properties of functions satisfying (A) and (B). In section 
4, we consider functions on products of topological groups, satisfying (A) and (B). 


2. PRELIMINARIES 


In this paper, the following notations and conventions are followed unless 
otherwise stated. 


G denotes a multiplicative group with identity e. A function on G means a 
complex-valued function. ‘Additive function’ on G is a group homomorphism of G to 
the additive group of complex numbers. If H is a subgroup of a topological group G 
then the topology on H is the subspace topology. C", R" have their usual inner pro- 
duct space structures. CG* denotes the multiplicative group of non-zero complex 
numbers. A non-zero function fon G means that fis not identically zero. We say that 
fonGisa D’Alembert’s function if it satisfies (A). For definitions we refer to 
Hewitt and Ross’. 


(This work was done when the author was a Teacher Fel!ow of the University Grants Commission 
at the Centre of Advanced Study in Mathematics, University of Bombay, Bombay 98.) 
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Remark 2.1 : A function f on G_ satifying (A) is non-zero if and only if 
f(e) =1. 

Remark 2.2: The set W = {z € C*: | z—1] < v3} contains only one subgroup 
of C* viz {1}. 

In Remarks 2.3, 2.4 and 2.5 below, let fon G be a non-zero function satisfying 
(A) and (B). In the proof of Theorem 2 in Kannappan‘, Remarks 2.3 and 2.4 below 
have been proved. 


Remark 2.3: Let f(x) = + 1 for all x € G. Then fis a homomorphism of G to 
the multiplicative group {—1, 1} and f(x) = (f(x) + [f(x)]~»)/2 for all x in G. 


Remark 2.4: Suppose there is Xp in G such that f(x.) ~A+1. Let B* = (/(x,))?—1. 
Then g defined by g(x) = f (x) + (f (xx0) —f (x) f (%o))/B is a homomorphism of G 
to C* with f(x) = (g(x) + [g(x)])/2 for all x in G. 


Remark 2.5: Let g = fif f(x) = + 1 forall x in G and let g be as in Remark 
2.4 otherwise. Then g is a homomorphism of G to G* such that f(x) = (g(x) + 
[g (x)]~*)/2 for all x in G. 


3. SOME PROPERTIES OF D’ALEMBERT’S FUNCTIONS 


Theorem 3.1—A function f on G is non-zero, bounded and satisfies (A) and (B) 
if and only if there is a character % of G such that f (x) = Re % (x) for all x in G. 


If G is a topological group, theorem holds when function, character are replaced 
by their respective continuous versions. 


PRooF : Suppose / satisfies the hypothesis, then g defined as in Remark 2.5is a 
homomorphism of G to G*. Since f is bounded 50 is g. Hence, | g(x) | = 1 for all 
x € G. Thus g is a character of G. Since 


F(x) = (8 (a) + [g ())/2 = (g (x) + 8 (%)/2, 
we see that f = Re g. Converse is obvious. 
If G is a topological group, f and g can be taken to be continuous. 


Remark 3.2: Ljubenova® has proved that O’Connor’s Theorem 5 holds for any 
abelian topological group. It follows from Remark 2.1 and above Theorem that 
O’Connor’s theorem holds when G is an abelian topological group and for non-abelian 


group if f satifies (B). Above argument shows that O’Connor’s theorem follows from 
Kannappan’s theorem. 


Proposition 3.3—Let H be 
group G. Let f on H be continu 
on G also Satisfying (A) on G. 


a subgroup of a locally campact Hausdorff abelian 
ous and satisfy (A), Then f has a continuous extension 
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Proor : If f= 0 on H, the zero map on G is the obvious extension. So let 
J = 0. Then by Corollary 1 of Theorem 2 of Kannappan‘, there is a continuous 
homomorphism g of H to C* such that f (x) = (g (x) + [g (x)]"’)/2 for all x in H. Let 
X (x) = g (x)/ | g (x) | and & (x) = log | g (x) | for x in H. Then % is a character and 
E is a real character of H and g (x) = exp (& (x)) % (x) forall x € H. Note that X and 
— are both uniformly continuous. Since 7 = {z € C:|z| = 1}, R are complete 
Hausdorff uniform spaces and H is dense in clH (closure of H), %,& can be uniquely 
extended to continuous functions %’, &’ over clH. It is easy to verify that %’ isa 
character and ©’ is a real character of clH. Further %’ and &’ admit continuous exten- 
sions %, and &, over G, where %, is a character and &, is a real character of G (see 
Hewitt and Ross’, pp. 380 and 391). Define g, and f; on G by g, (x) = exp (8 (x))%: (x) 
and f, (x) = (g; (x) + [g; (x)]"’)/2. Clearly f, is the required extension of /f. 


Remark 3.4: Let H and G beas in Proposition 3.3. Let L be a continuous 
additive function on H. Then L has a continuous extension to G, which is an additive 
function on G. 


Proor: Apply the same technique as above extending L first to clH and then 
to G. 


Proposition 3.5—The following are equivalent on G: 


(i) For a non-zero function f on G satisfying (A) and (B), there exists an additive 
function L on G such that f (x) = cos L (x) for all x in G. 


(ii) For any homomorphism g of G to C*, there is an additive function L on G 
such that g (x) = exp (iL (x)) for all x © G. 


(iii) For any character % of G, there is a real character — of G such that 
7%. (x) = exp (i & (x)) for all x in G. 


If G is a topological group, Proposition holds when function, homomorphism, 
character, real character are replaced by their respective continuous versions. 


Proor: (i) > (ii): Let g bea homomorphism of G to C*. Then f (x) = (g (x) 
+ [g (x)}-)/2, x € G satisfies (A) and (B). So by (i), there is an additive function H 
on G such that f (x) = cos H (x) = (exp (iH (x)) + exp (—iH (»)))/2 for x € G. 
Hence, by Theorem 3 of Kannappan’, g (x) = exp (/H (x)) or g (x) = exp (— iH (x) 
for all x in G. So take L = H or —H. 


(ii) > (iii): Given a character % of G, by (ii) there is an additive function 2 on 
G such that % (x) = exp (iL (x)) for all x in G. Clearly L is a real character since 
Im L (x) = 0, | % (x) | being | for all x. 

(iii) > (i): Suppose / satisfies (A) and (B) and is non-zero. By Theorem 2 of 
hism g of G to C* such that f(x) = (8g (x) + 


an’, there is a a homomorp 
eo) = g(x)! | g (x) | and ¢ (x) = log | 8 @) | for 


[g (x))°')/2 for all x in G. Let % (x) 
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x € G. Then & is a character, ¢ is a real character of G and g (x) = exp (¢ (x)) % (x). 
By (iii), there is a real character & of G such that % (x) = exp (i & (x)). So g (x)=exp 
(¢ (x) + i & (x)). Let L (x) = —i (¢ (x) + 7 & (x)), x EG. Clearly L is as stated in (i). 


If G is a topological group, each of functions involved above can be taken to be 
continuous. 


Theorem 3.6—Let f be a non-zero, real-valued function on G satisfying (A) and 
(B). Then : 


(i) fis bounded if and only if f(x) = Re % (x) for some character % of G. 


(ii) f is unbounded if and only if f(x) = ¢ (x) cosh — (x) where ¢ is a homo- 
morphism of G to the multiplicative group {1, — 1} and = is a non-trivial real 
character of G. 


If G is a topological group, theorem holds when function, homomorphism, 
character, real character are replaced by their respective continuous versions. 


PRroor : (i) is Theorem 3.1. 


(ii) f (x) = (g (x) + [g (x)])/2, where is a homomorphism of G to G*. Since 


f(x) is real, for every x, | g (x) | = lorg (x) is real. So, range of g is contained in 
X = {z:|z| = 1} U {non-zero reals}. Since g isa homomorphism, range of g is a 
subgroup of C* contained in Y. But then it has to be either contained in {z:] z] = 1} 


or {non-zero reals}. 


If g (G) C {z:|z| = 1} then g is character of G and hence for any x E G, 
J (x) = Reg (x). So f is bounded, which is a contradiction tof is unbounded. Hence 
g(G) must be contained in {non-zero reals} and | g(x) | # 1 for at least one x E G. 
So g is a non-trivial homomorphism of G to the multiplicative group of non-zero 
reals. Let ¢ (x) = g (x)/ | g (x) | and & (x) = log | g(x) | forx EG. Then ¢ is a homo- 
morphism of G to {—1, ]}, E is a non-trivial real character of G and g (x) = exp 
(E (x)) ¢ (x). Hence for any x € G, 


J (x) = (exp (& (x) $ (x) + [exp (E (x)) ¢ (x)]})/2 
= # (x) (exp (E (x)) + [exp (E (x))}-")/2 = ¢$ (x) cosh & (x). 


The converse is obvious. 


If G is a topological Stroup and f is continuous, then by Corollary 1 of Theorem 
2 of Kannappan‘ there is a continuous homomorphism & Of G to C* such that 


f(x) = (8 (x) + [g (x)] "D2 for all x EG. H = 
& (x) = log | g (x) | are continuous on G, Ravers C18 Cie (iran 


Remark 3.7 : If G is such that G2 == G and if ¢ j ; 
Me. = nd if ¢ is a homo 
the multiplicative group {1, iy thee : morphism from G to 
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D’ALEMBERT’S FUNCTION ON PRODUCTS OF TOPOLOGICAL GROUPS 


In what follows, if {G,: € J} isa family of (topological) groups then M1 G, 
denotes the direct product (with the product topology) and I1* G, denotes the weak 
direct product (with the topology induced as a subspace of the direct product [I G,; 
with the box topology) unless otherwise stated. For a summation or product that 
follows in this section index set is J unless otherwise stated. 


We identify each G, with the subgroup of II G, formed by the x’s all of whose 
coordinates except the jth are equal to the identity. Thus for each / € J, G; C II G, 
and each x, is an element of II G;. And similarly for Il* G;. 


Theorem 4,1—Let G = II G,, where each G,,/ € J is a topological group. Let 
f on G be non-zero and satisfy (A) and (B). Then the following are equivalent : 


(i) f is continuous. 


(ii) f (x) = (II g, (x)) + [Il g, (x)}})/2 for all x (= < x; >) in G, where 
{g,: G,; > G*, 7 € J} is a family of continuous homomorphisms in which only finitely 
many homomorphisms are non-trivial. 


Proor : (i) > (ii): By Corollary 1 of Theorem 2 of Kannappan*, there is a 
continuous homomorphism g of G to G* such that f(x) = (g (x) + [g (x)]))/2 for 
all x in G. Let g, be the restriction of g to G,. The g),/ € Jisa continuous homo- 
morphism of G, to C*. Hence, assertions (1) and (2) given below are sufficient to 
prove the result. 


(1) g; is identically equal to 1 for all but finite number of /’s and 
(2) g (x) = I g, (x,) for all x in G. 

Let W = {z € G* :|z — 1| < v 3}. Since g : G > C* is continuous, 
U = g~ (W) is an open set in G. Hence by definition of product topology, all but 
finite number G,’s lie in U. Consequently, for all but finite number of j’s, g (G;) C W. 
Hence using Remark 2.2, we have, for all but finite number of j’s, g(G,) = {I}. 
Accordingly g; = g | G; =! for all but finite number of j’s. Assertion (2) can be 
easily checked. 

(ii) > (i): Let {g,: G; > C*,/ € J} and f be as in (ii). Suppose 8, is identi- 
cally equal to 1 for all j’s except /1, Jz ..--- _ ine Define g : G -> C* by g (x) = Ig /(x;), 
x= <x, >€G, Since ge (x) = gyn (xj), K = 1,..5 7 are continuous and 


» oe . . 
g (x) = Il & (x), g is continuous. Hence f is continuous. 
k=l 


Definition 4.2—We say that a topological group G satisfies condition (C) if for 


i i i G such that 
every continuous character % of G there is a continuous real character E of 


% (x) =exp (i & (*)) for all x in G. 
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Remark 43: If His a subgroup of R", then H satisfies condition (C). If each 
G, in Theorem 4.1 satisfies condition (C), then g; (x,;) = exp (i L; (x,)) where L, is 
an additive homomorphism of G; to C. So we have Theorem 4.4. 


Theorem 4.4— Let G = I1 G;, where each G,, j € J isa topological group 
satisfying (C). Let fon G be non-zero and satisfy (A) and (B). Then the following are 


equivalent : 
(i) f is continuous. 


(ii) f (x) = cos [2 LZ, (x;)] for all x in G, where {L, : G,;>C,jE€ Jsisa 
family of continuous additive functions in which only finitely many are non-zero. 


(iii) f(x) = cos L (x), x € G for some continuous additive function Z on G. 


Lemma 4.5—Let H be a subgroup of the additive groug R" (with the usual 
topology). An additive function Z on H is continuous if and only if there is a vector 
a in C” such that Z (x) = a. x for all x in H, where dot denotes the usual inner pro- 
duct in G”. If L = 0 then we can take a = 0. 


Proor: If Lis continuous, by Remark 3.4, L has a continuous extension L’ 
over R" that is a continuous linear function of R" to C. Hence, there is a vector a in 


C” such that L’ (x) = a.x for all x in R". Consequently, ZL (x) = a.x for all xin H. 
Converse is obvious. 


Remark 4.6: Let H be a subgroup of R". Then using the technique of the above 
lemma it is easy to see that a real character — of H is continuous if and only if there 
is a vector ain R*” such that & (x) = a.x for all x in H, where dot denotes the inner 
product in R". 


As an immediate consequence of Remark 4.3, Theorem 4.4 and Lemma 4.5 we 
have the main result of this paper. 


Theorem 4.7— Let G = II G;, where each G,, 7 © Jisa subgroup of R*" (with 
the usual topology). Let f on G be non-zero and satisfy (A). Then f is continuous if 
and only if f(x) = cos (= a;.X,) for all x in G, where (aj: / € J} isa set of vectors in 
CG”, all but finitely many of which are zero and dot denotes the inner product in C’, 

: Remark 4.8 : Let {G,:7€ J}bea family of abelian topological groups. Let 
Il* G, be the weak direct product with the subspace topology of I G, with the pro- 
duct topology. Then TI* G, is a dense subgroup of Il G,;. Hence, using technique 
of Proposition 3.3 and Theorem 3 of Kannappan‘, it can be easily shown that a con- 
sane function f on ue G, satisfying (A) has a unique continuous extension f* to 

G, which satisfies (A). Using this it is Casy to see that Theorem 4.7 is valid if we 


ae G as the weak direct product n* G, of subgroups G, of R" with the topology 
induced as a subspace of 1G, with the Product topology. 
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Now we consider functions satisfying (A) and (B) on the weak direct product of 


groups. We state the following result, which can be easily proved by using Theorem 
2 of Kannappan‘. 


Theorem 4.9—Let G = II* G,;, where each G,, 7 © J is a group. Then f on G 
is non-zero and satisfies (A) and (B) if and only if there is a family {g,: G > C*, 
J € J}of homomorphisms such that f(x) = (Il g, (x,) + [II g; (x,)]7?/2 for all x 
in G. 


Recall that, from now onwards, if {G,:j € J} is a family of topological groups 


then II* G, denotes the weak direct product with the topology induced as a subspace 
of the direct product II G; with the ‘box-topology’. 


Theorem 4.10—Let J be a countable index set and let G = I* G, where each 
G,;,j © J isa topological group. Let f on G be a non-zero function satisfying (A) and 
(B). Then the following are equivalent : 


(i) f is continuous. 


(ii) There is a family {g, : G; > C*, j © J} of continuous homomorphisms such 
that f(x) = (II g, (x,;) + [II g, (x,)]-)/2 for all x in G. 


(iii) There is a continuous homomorphism g of G to C* such that f (x)=(g(x) + 
[g (x)]~1)/2 for all x in G. 


ProoF : (i) > (ii): By Corollary | of Theorem 2 of Kannappan‘, there is a 
continuous homomorphism g of G to C* such that f (x) = (g (x) + [g (x)]~’)/2 for all 
x € G. Since g is a continuous homomorphism of G to C*, the restriction map g, of 
gtoG,, 7 € Jisacontinuous homomorphism of G, to C* and g (x) = II g, (x,) for 
all x € G. Hence, {g; = g| G,:j € J} is the required family. 


(ii) > (iii) : Let {g,;: G; > C*, 7 € J} and f be as in (ii). Then g(x) = II g, (x,), 
x € Gisa homomorphism of G to C*. Hence, it only remains to show that g is 
continuous. 

If J is finite, then g is obviously continuous. So, let J = N. Lete > 0 be given. 
Since g, is continuous, U; = ieee] By Ly Li «/2} is a neighbourhood of e in 
G,. Let U = (11 U,) - (II* G)). Then U is a neighbourhood of the identity in the 
subspace topology of G induced as a subspace of II G, with the box topology. If x is 
a point in U then for any N € N, 


N N-1 Nw 
|Mg,(x)Jy-—I!1s | o g, (x) || gn (%w) -— 1/1 4+ | T 8; (x,) | 
1 1 
| Zv—s (Xy-1) | # oe + | 81 (%1) — II 


on cer oF ag 
< fia + 2B le) — 11) < et Ele @) 11. 
1 1 
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Hence for any x € U, 
Co 
CE be =|Megyx)—-1|<e2laO)~loes 


So, g is continuous at the identity and hence continuous on G. 
(iii) = (i) is obvious. 


Proposition 4.11—Let / be a countable index set and let G = II* G;, where each 
G;, € Jis a topological group satisfying (C). Let fon G be a non-zero function 
satisfying (A) and (B). Then the following are equivalent : 


(i) f is continuous. 


(ii) There is a family {L):G, > C,j € J} of continuous additive functions 
such that f (x) = cos [2 L, (x,)] for all x € G. 


(iii) There is a continuous additive function L on G such that f (x) = cos L(x) 
for all x in G. 


Proor : Using Theorem 4.10 and Proposition 3.5 the result follows. Note that 
the sum &Z,(x,) has only finitely many non-zero terms as any element x in G has 
only finitely many coordinates different from identity. 


Corollary 4.12—Let {G,:j € J} be a countable family of subgroups of R" (with 


the usual topology) and G = II* G,. Let f on G be non-zero and satisfy (A). Then the 
following are equivalent : 


(i) f is continuous. 


(ii) f (x) = cos (2 a;.x,) for all x in G, where {a,:j € J} is a set of vectors in 
C” and dot denotes the inner product in C”. 


Remark 4.13: In Theorems 4.1, 4.4, 4.7, 4.9, 4.10, Proposition 4.11 and 
Corollary 4.12, if f is also bounded then using Theorem 3.1 and methods used in 
respective results it is easy to show the following : 


(i) In Theorems 4.1 and 4.10, g,'S are continuous characters and f(x) = Re (IIg; 


(x,)). 
(ii) In Theorem 4.9, g;’s are characters and f(x) = Re (II g, (x;)). 


(iii) In Theorem 4.4 and Proposition 4.11, Z,, 7 € Jand L are continuous real 
characters of G,; and G. 


. (iv) In Theorem 4.7 and Corollary 4.12, a,;’s are vectors in R" and dot denotes 
inner product in R®. 


Remark 4.14: In Theorem 4.10, Proposition 4.11 and Corollary 4.12 if the 


index set J is uncountable then (i) > (ii) but (ii) = : se 
following : (ii) =/> (i). This is shown by the 
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| Counter example ; Let J be uncountable and let G = I* G,;, where each G),=R 
(with the usual topology). Let &’: R — R be defined by &' (t) = at for allt E R. 
For each j € J, let §&; = &’ and %; = exp (i &’). Then {€,: 7 € J}isa family of con- 
tinuous real characters and {%,: © J}is a family of continuous characters. Define 
&, Zon G by & (x) = DE, (x,) and X(x) = mm X,(x,;) = exp(id E, (x;)). 


Claim : =, % are discontinuous. 


Let U = (Il U,) ( (I1* G,) be a basic neighbourhood for the identity in the subspace 
topology of G induced as a subspace of II G, with the box topology. Then {U pee Aes 
is an uncountable family of neighbourhoods of zero. So there is a natural number my 
such that (—1/no,1/n,)CU, for each jin some uncountable subset A of J. Let a,, Gs; 
Gno4, be any yo + 1 elements in A. 


eeey 


Let x = < x; > inG be such that 
xj = 1/n, + l if 7 = a), Aayeeey Qn0+1 
= 0 otherwise. 


Then x € Uand & (x) = 7 and X (x) = —1. Thus any basic neighbourhood of the 
identity in G contains x such that § (x) = 7 and % (x) = —1. Hence &,% are discon- 
tinuous. Since ¥% is discontinuous, f = Re X is also discontinuous. Thus (ii) =/=> (i). 


Remark 4,15 : Kannappan® has shown that his condition 
S (xyz) = f (xzy) for all x, y, z in G (B) 


is very restrictive and reduces the study of (A) on a non-abelian group G to the study 
of (A) over the abelian group G/G’, where G’ is the commutator subgroup of G. In 
view of this it is sufficient to study (A) on an abelian group. However, we study 
functions satisfying (A) and (B) on non-abelian groups because the proofs of the 
results proved in this paper are same in both the cases. Also we have advantage of 
studying (A) on the original group instead of studying it on the quotient group. 
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The primary purpose of this paper is to show how techniques and topics in 
applicable analysis are put together to obtain the result. 


There is a plethora of literature on the gamma function; hence, an explanation for 
this paper is in order. The result, lim I (iy) as y ~ co = 0 is not the motivating 
factor for writing this paper although the result has an application for the calculation 
of the absolute value of factorials of complex numbers. The main result was given by 
Stieltjes’ in 1889 but the technique used is not accessible to the average reader, Thus, 
the primary purpose of this paper is a nice form of exposition to show how techniques 
and topics in applicable analysis are put together to obtain the result. 


Techniques and topics used (quite different than in the work of Stieltjes) follow 
this sequence : (a) the extraction of a factor from an infinite product, (b) the repre- 
sentation of an infinite product by two infinite products, (c) the examination of conver- 
gence of an infinite product by determining an upper bound for a sequence of partial 
products, (d) the use of Euler’s product fo nz/sin 7z, (e) the use of hyperbolic func- 
tions and (f) the computation of | z! |. These techniques and topics should be in 
the repertory of the applied analyst. 


In some texts?’? the following problem is given as an exercise: 
Prove that 


| C (yi) | *? = z/y sinh ay. 
The student is expected to use the reflection formula for the gamma function 


7 
— Z = ————_—— 
B (2) eS ) sin 7Z 
However, the reflection formula. It is only fair to state, is not a simple matter to 
obtain. A detailed development of the above formula can be found in Ross*. Our 
starting point is Weiersirass’s formula for the reciprocal of the gamma function 


io) 


re mse [](1+5)°" ons 


a | 
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where y is Euler’s constant. Let z= x + yi. We note that 





| 2 2 1/2 
slojie de Aa ((r4 sy 48) 
jit 2] art Se th n® 
and that | eY*+ | = eY*, Egn. (1) can then be written as 
+ 2 2 \1/2 
l = (x= LC y2\t/2 eY* I (( l os x ) + ea ) esl”, 
Rice yew ee 1h . : 


...(2) 


The factor ( nes is taken out of the infinite product above. Motivation 
n 


for this step is explained later. The infinite product in eqn. (2) can be represented by the 
product of two infinite products : 


ee | ee E 


n=1 


The first infinite product above is the infinite product in Weierstrass’s definition 
eqn. (1) for 1/I’ (x) and its convergence is well-known. Of immediate concern is the 
convergence of the second infinite product. One can, of course, write out a few terms 
of the infinite product in eqn. (2) and a few terms Of the products in eqn. (3). 


From the identity that ensues, one cau intuitively assume convergence of the 
second infinite product. However, we shall not rely on intuition; our recourse will be 
on the reliance on theorems of convergence of sequences. For convenience in writing, 
let y*/(x + n)* be denoted by w,. 


We observe that for fixed x and y, the sequence of partial products (1 + u,)'!2, 
(1 + uw)? (1 + uw), (1 + wy)? a + uy)? (1 + us )'2, .. isa sequence of positive 
monotonically increasing functions of n and, accordingly approaches a limit if bounded. 
To determine if the sequence is bounded we write 


u 
l+u,sSe" 


so that 


2 
. 


(l+u)ji2sen 
Then 


CL + ay)" (1 + uy)? (l +4,)!? < el Pe ie “al” 


Now let » tend to infinity and we will have 
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coo co 
uy (1 + um)'/? S exp (4 = u,). 
bdo ] 


oo 
But a u, is convergent. Thus, the infinite product has an upper bound. Because the 


sequence of partial products is monotonic increasing and has an upper bound, the 
infinite product is convergent. 


The next step is to put the infinite products in eqn. (3) in closed form. The first 
infinite product is the same product that appears in Weierstrass’s definition in eqn. (1) 


so that 
co 
sid —tin l 
Hti+ =): coh In bt 5 
n=) 


When the expression on the right above is substituted into eqn. (2) we will have 
l ie | 2 \IR 
Tre +o teso I (1+ eta) - te 
The squaring of both sides and taking reciprocals gives 
00 
res pi>> SEEM TT (1+ Far) - ...(5) 


x? + y? (x-- 7)" 


n= 


An equivalent and simpler form of the above can be obtained if we start the 
product with n = 0: 


iPro t+ yi) |? = (FQ) Il ( ie cor) - 6 


n=0 


Physicists use the above relation in calculations in the theory of beta ray decay. 
To evaluate I’ (yi), let x = 0 in eqn. (5) and insert i? = — 1 in the denominator of the 


infinite product getting 
co ye ; 
ae — (yi 7 
Ironi" = “all ( ea ( 
imi 


The purpose of extracting the factor ( | . ) in eqn, (2) is now clear. The 
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infinite product in eqn. (7) can be put into closed form. We recall Euler’s infinite 


552 


tN 
02 





nr 


product for sin 7z is 
) 


=~ M(-: 
...(8) 


sin 7Z 
n=) 








With z = yi, the product above is 
myi Ty Bs my 
sinzyi  —s 7 sinh zy sinh my 
When the above is put into eqn. (7) we have the result 
SAC). 


7 
(yi) | = ————_—., 

|T Qs) | Tinka 
If we treat co as a number, we can write 


[' (co i) = 0, 
The result (9) is the same for argument —yi, as can be seen in (8), because 
= 2 + 4i is to use Legendre’s duplication 


sinh (—nyi) = —sinh zyi. 
One method to calculate | z! | where z 
formula for I (z) getting 
(3. + 4i) = 0.005225... — 0.17254... 7. 
From elementary proper- 


Then 


However, it is more expendient to use the result is eqn. (9). 
1(2+ 4)! | =| P38 + 4] = [24 4/1 + 424i P(4i)]. 
. We will have then 


[(2 + 4) | =1 FG +41) = 0.1726... 
ties of the gamma function 

Equation (9) with y = 4 gives | T (41) = 0.00234049 ... 
| (2 + 4i)! | = (4/20) (v17) (4) (0.00234049) = 0.1726 ... . 


A useful by-product of this analysis stems from eqn. (6) by first writing it as 


co 
ey eae s. ) ip 
II ( Lar IEC) 
we can write 


nel 


[P+ yi)|* _ 
| P(x) |? he 


For fixed x and y the infinite product is less than unity. Thus, 


IT +y)1 S| PQ]. 


—s 
. 
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A symmetrizer of an abritrary matrix A is a matrix X satisfying the matrix 
equations XA = At¥ and X¥ = X'. A symmetrizer is useful in transforming 
an apparently nonsymmetric eigenvalue problem into a symmetric one. An 
algorithm to compute a matrix symmetrizer is presented and illustrated. 


1. INTRODUCTION 
A symmetrizer of an arbitrary square matrix A isa symmetric solution X that 
satisfies XA = A‘X, where f¢ indicates the transpose. Symmetrizers reduce a non- 
symmetric (real) eigenvalue problem (nep) into a symmetric eigenvalue problem (sep) 
which is relatively easy to solve and are useful in the stability problems of control 
theory’. These are also useful in the study of general matrices®. 


Transformation of nep to sep using a symmetrizer—Let A be a square nonsymme- 
tric real matrix and X be a nonsingular symmetrizer of A. Then the nep Ay = Ay can 
be reduced to an sep Bx = Ax as follows. Ay = Ayimplies YAy = AXy. Since X is 
real symmetric, we have X = P! DP, where P is orthogonal (i.e., P’ = P~) and 


D = (dy) is diagonal. Therefore, we have PYAy = ADPy, Let D, = (di, ) be a dia- 
gonal matrix such that Dj = D. Then PYAy =AD? Py. di + O for all i since ¥ has 
nonzero eigenvalues (as X is nonsingular). Hence D,* exists. Thus, Bx = Ax, where 


the symmetric matrix, real or complex, 


=1 t =i 
B = D;' PX AP' D, At) 


and x = D, Py. If X is positive definite then dj; > 0 and consequently dj, is real. In 


this case B is real symmetric. If Y has a negative eigenvalue then B is complex 
symmetric. 


Symmetrizer to compute zeros of a real polynomial—A symmetrizer can be used 
to obtain the roots of a polynomial equation® 


The associated companion matrix® of 
the polynomial 


* 


supported by CSIR. 
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Pr (x) — = qj x? 
—s 


where a, = 1], and a, i = 0, l,...,n — l are real, is 


fO 0 eee 0 =o 7} 
l = 
ita | Fos. & a | 
: | 
LO 0 tee l Bens | 


whose eigenvalues are the zeros of p, (x). To obtain these eigenvalues, we compute 
an X such that YA = A‘X and X¥ = X‘ and then we compute the eigenvalues of the 


symmetric matrix B = D,’ PXAP' D,', by a method, say, the Jacobi method*”” if 


the matrix X is positive definite. If the nonsingular X¥ has one or more negative 
eigenvalues then a Cholesky-type decomposition®’'? can be used to obtain the eigen- 
values of B. 


Existence of a symmetrizer—There exists a nonsingular symmetrizer for any 
square matrix’®. Let A be areal nonsymmetric matrix and X be one of its symmetri- 
zers. Then the following table depicts the scope of the symmetrizer YX. 


TABLE I 


Scope of a symmetrizer 


EE  ————— 


Eigenvalues of X Reduction of nep to sep Utility 
EN a a os ee ee ee ee eee 
All positive (X is positive Possible Useful since sep is easier to 
definite) solve (say, by Jacobi method) 

One or more zero Not — 
(X is singular) possible 
One or more negative Possibe Complex arithmetic needs to 


be used to solve sep (say, by 
Cholesky-type decomposition). 


SO 


Since X is not unique, one can always compute another X by a different choice 
c. 2) so that X is positive definite, a positive 


and no zero 


of some of the elements of X (see se 
definite ¥ does not, however, exist if the matrix A has complex eigenvalues. 
2. ALGORITHM TO ComPUTE A MATRIX SYMMETRIZER 
Let A = (aij) be ann xn real nonsymmetric matrix. Also, let X =(x1;) bea 
symmetrizer of A to be computed. 
Step 1: (Compute Au) 


Fori: = l1ton — 1 do 
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compute the (n — i) X (n — i + 1) matrix Aji, 
where (k, t) th element of Ay 
[= Gui kt anift +i —-1=k + ielse 
> = e¢i-1,k4t- 
Step 2: [Compute Ai, j > i] 
Fori: = 1 ton — 1do 
compute the (n — i) x (mn — j + 1) matrix Aj, 
where (k, t) th element of Aj; 
:= O0fork:=Iltoj—i-1 
>= Ofork:=j—i+lton—iandt4#k -—j+i4+1 
>= — Ayyt-1y fork: =j —i 
>= —ayfork:=j—i+1ton—iandt=k—j+i+1. 
Step 3: [Compute Aj, j < i] 
Fori: = 2tom — 1do 
compute the (n—i) x (mn — j -+ 1) matrix Aj), 
where (k, 1) th elememet of Aj, 
:-= Ofort: = 1 toi—j 
>= Ofort:=i-—j+2ton—j+ landk+~t+j-—i-1 
>= 444 fort: =i-—j+1 
= — an fort: = i—j+2ton —f+ 1 andk 2+ j7— }— 1, 
Step 4: [Obtain C] . 
Obtain C = (4),)i = 1, 2,...,n — Pe a hy Soy 
Remark : C is an (n? — n)/2 x (n® + n)/2 matrix. 
Step 5: [Solve Cz = 0] 
Solve Cz = 0, 


where 


Z_= [x11 X10 eee Xin Xo» X53 eae Xon ees Xn=t9n-1 Xn-19n eel 


as follows. 


SYMMETRIZING A MATRIX 557 


Step Sa: Choose, if none of the columns of C is zero, m elements x;; = c + 0, 
x12 = *13 =... = Xi, = 0 or any n elements x;; (i = j) so that the resulting equations 
are consistent. Ifa column, which constitutes the coefficients of some x1), of C is 
zero then choose x;, along with any other n — | elements of X so that the system is 
consistent and at least one of these n — 1 elements is nonzero to ensure nonsingularity 
of X. Ifk columns of C are zero then choose the corresponding x1, along with any 
other n — k elements of X (without disturbing consistency) so that at least one of 


these n — k elements is nonzero to ensure nonsingularity of Y. Let the resulting equa- 
tions be C’ z’ = b’. 


Step 5b: Use any method, say, the Gauss elimination method*” to solve the 
equation — this gives Y and terminate. 


Proof of the algorithm—The proof steps 1 — 4 follows from the expansion and 
rearrangement of the equations YA = A'X, X' = X and that of Step 5 (solving linear 
equations) is well-known. 


2.1. Stability of the algorithm—The algorithm comprises (a) computing C, which 
is always stable, and (b) solving the linear equations C’ z’ = b’. Therefore, the 
numerical stability of the algorithm essentially depends on the degree of singularity of 
the matrix C’ and on the stability of the method that is employed to solve C’ z’ = D’. 
The choice of the elements in Step 5 decides whether C’ is near-singular, i.e., ill- 
conditioned with respect to the inverse, or not. There is, however, no easy way to 
know beforehand which will be a good choice. The selection of a stable method to 
solve linear equations usually poses no problem. 


If the matrix whose symmetrizer ¥ = (x1;), where X,, Xp,...,X, are the rows of X, 
is required is lower Hessenberg and if we choose x, of the form[c00... 0], where 
c ~ 0 (see Sec. 4 also) then the algorithm is stable provided the Hessenberg matrix has 
the codiagonal nonzero. If any one element of the codiagonal is near-zero then the 
algorithm tends to become unstable unless some other appropriate choice is made. 


3. EXAMPLE 


We compute a symmetrizer of the matrix A = (aj,), where 














ym me foe 
PPA) a a a 
1-1 2 
Here 
1 —-2 —1/|0 -! 0 
Cwi 0 0-3/0 0-1 | 
J 
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and 
Z v= [X11 X12 X18 Xoo Xo9 X33)". 
Choice 1—X2, = 3 (this choice is made since the corresponding column is null), 
X;, = 2, and x,, = 1/2. Any other choice for any three elements of X may be made 


preserving the consistency of the equations. Choice | gives the nonhomogeneous 
equations C’ z’ = b’, where 


—|] —-] 0 X13 - —] 
c=] —-3 O-1 [2 =| 3 |b’ = 0 
—l —1 ] X33 uU Qu. 

Applying the Gauss reduction method with partial pivoting and then using the 
back substitution, we obtain x33 = 1, X,3 = 4/3, x,, = —1/3. Hence the symmetrizer 
is 

2 1/2 —-1/3 
om 1/2 3 4/3 
—1/3 4/3 l 


which is positive definite. The symmetrizer transforms nep Ay = Ay to sep Bx = Ax, 
where B = D,' PXAP! D>’, X = P' DP, and D = D? for some diagonal matrix D,. 
From X, we compute using the Jacobi method 


. 1787430 8919905 .4152155 
P =| —.9434132 0355713 .3297064 
.2793252 —.4506525 .8478737 


7 3.720851 0 0 
D= 0 2.097642 0 
0 0 0.1815073 


Taking the positive square-root of the diagonal elements of D, we obtain D, (negative 
square-root can also be taken). Hence 
2.7547790 —1.1969110 —0,1342376 
B=] -—1.1969110 4.3387318 —0.4728604 
—0.1342376 —0.4728604 2.9302057 


It can be easily checked that the eigenvalues (computed by the Jacobi method) of B, 


viz., 5.042097, 2.999203, 1.982418 are approximately the same as those of the real 
nonsymmetric matrix A, which are = ee nef 


Choice 2—X2s = 1,%1 = 1, x. = 0, Solving the resulting nonhomogeneous 
equations, we obtain the nonsingular symmetrizer 
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1 0 —1/3 
X= | as 4/3 
—1/3 4/3 l 


having a negative eigenvalue, viz., —0.3743685. Hence B = D,* PXAP! D,* 


1.8348370 —0.2159493 10.4159453 
= —0.2159492 5.1176470 i0.5438464 
i0.4159453 = 10.5438464 3.0475158 


where 
2.3743690 0 0 
D= 0 1 0 . 
0 0 —0.3743685 
and 


0.1714986 —0.6859943 —0.7071068 
P =| 90.9701425 0.2425356 —0.7304784E-20 
0.1714986 —0.6859943 0.7071068 


The eignvalues of the complex symmetric matrix B can be computed by a Cholesky- 
type decomposition; however, it can be easily checked that the eigenvalues of B are 
approximately the same as those of A, viz., 2, 3, 5. 


All the foregoing computations have been carried out on DEC 1090-TOPS 10 
computing system in single precision and 7 decimal digits are input. 


4. Recursive ALGORITHM FOR A LOWER HESSENBERG MATRIX 
A simple and elegant recursive algorithm to compute a symmetrizer was sug- 
gested by Datta®. The algorithm is as follows. Let A = (a,;) be ann x n lower 
Hessenberg matrix with nonzero codiagonal. Let 1, X2,...5 %n be the rows of a 
symmetrizer X to be computed. 


Step 1: Choose x, # 9 arbitrarily. (To ensure nonsingularity of X, choose x, 
of the form [c 0 0...0] where c = 0.) 


Step 2: Compute X,-1, %n—29- ™1 recursively from 
” Sse (1/4114) (X14: Aim Alpisty, X41 — Alsosl41 Xie — --- Ansley 3 Fh 


This algorithm Is applicable only to lower Hessenberg matrices with nonzero codia- 
gonal elements since each codiagonal element 4,i+1 OCCUrS in the denominator of 
Step 2 of the algorithm. However, any matrix can be reduced to a lower Hessenberg 
form by using similarity transformations involving Gauss-type elimination and this 
reduction is direct; thus, the eigenvalues remain invariant. 
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Derivation of the recursive algorithmn—The foregoing recursive algorithm is a 
special case of the suggested general algorithm and is derived as follows. Consider 
the homogeneous equations Cz = 0 for then X mnonsymetric real matrix 4d = (414), 
where C and z are defined in Sec. 2. Substitute a,; = 0 for 7 = i + 2,7 + 3,..., nin 
C since A is assumed to be lower Hessenberg. Choose (as suggested by Datta) 
Xuq = Xp = C FO, Xyg = Xen =O, ..-5 Nee = O in Zz and) write the Keauling linear 
equations C’ z’ = b’, where C’ isa triangular matrix. Using back substitution, the 


recursive algorithm follows. 


If A has any other special property then this can be incorporated in the general 
algorithm so that the general algorithm reduces to a simpler (special) one. 


5. CONCLUSIONS 


A general algorithm for computing a symmetrizer of an arbitray matrix A 
is suggested. A case where A is lower Hessenberg is shown to be a special case of the 
general algorithm. Ifthe nonsymmetric matrix has any other special structure then 
this structure can be exploited to generate aspecia! case of the algorithm. The 
advantage of obtaining a symmetrizer for computing the eigenvalues of a nonsymme- 
trix is yet to be fully investigated. 


The computation of a symmetrizer involves a finite number of arithmetic opera- 
tions and hence a symmetrizer can be computed exactly using finite-field transform 
techniques*"*’*’? superimposed on the (direct) method of solving linear equations. The 
need to have an exact symmetrizer arises when errors involved due to the use of real 
arithmetic is considerable. 
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In this paper, certain symbolic relations and expansions associated with basic 
hypergeometric functions of three variables have been established. 
1. INTRODUCTION 


Following Burchnall and Chaundy’”’, Srivastava‘ has obtained single and double 
series expansions associated with hypergeometric functions of three variables by 
employing the inverse pair of operators : 


f h, h - é 3 8, + 33 } =. aa _ > 
Vxsysz (h) — E: | | ao DS ( Sr+s ( 55)r ( S3)s 
J 





(A+ 5,A4+ 5, + 3; 30 (A)r+s Ur ()s 
f h+o6, A+ do + 85 1 
Ax,y5s (/) th | | 
Lh, h+6,4+5,4+ 83 j 
Si a wed) 27 Gok hed oe Q2) 
(1 —h—-3$,— os — 8 ras l r (1 s - ) 
r,sP0 . ( ) ( ) 
where 5 x od 5 2 , ° 
C — Y =—- Ca = | are ae } ae 
1 ox , 2 J oy > 3 “ az " 


Here we shall establish certain symbolic relations and expansions of basic hyper- 
geometric functions of three variables by making use-of the following g-analogues of 
the inverse pair of symbolic operators (1) and (2): 


(A+, } 
vi (A) =m | yer Stes (3) 
a Lhh+o+¢e+u | y 
Cs) 
A hye a | rh+o+get+y is i 
Xyez 1 h + 6, h 4. d a v J ei 


The series equivalents of the Operators in (3) and (4) are given by 


and 


where 


and 
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Vere @) > (9): (Q)s r+ ses) 
r,sP0 
AW@ = (—8)u; (—¢), (—¥), gore 
ee 2 @ @). -h-0-$-Dre, 6) 
Mek 
é é 0 
a Sa any By ,u = 2 , g® = exp (6 log q) pany} 


(An = (1 — 6) (1 — 6q) ... (1 — 6q"-") ; (6)o = 1 


The basic analogue of Srivastava’s hypergeometric function of three variables® (p. 428) 
is defined as (see also Srivastava and Karlsson®, p. 44) : 


(Ay a2(6)5 (o)3'():2 (8) (7).5'(8) : } 
q,;X,),z | 
18 dese ay al C8 Fe Ct Beat Col De OF hd TA a ' 
[(a) mint? [OV min (On+2 (dram [(e) Im [CF Jn 
[(2') m+n+p [(b')Jntn [(c’)]nap [(4')]o4m [(e’)] (Cf a)P 


m,n,pPO 


pe [(g)lo x™ y” 2? 
[(g’)lo (9)m (Qn (Qo 


(3) 


co eee 


= 


...(8) 


2. THE SyMBOLIC RELATIONS 


Employing the operators defined by (3) and (4), we have the following symbolic 


relations : 


F as: —;—;—:),; ba; bs; ! 
os) | Q;X,y.2 | 
Up etl ett Beek Pet Os J 
= V2, (2, (2, by; c15 x). OM (G, ba, B35 Cay C35 Ys 2): -..(8) 
f a::—3;—; — 2,3 by; bs; } 
Oo | Gs%,¥,Z | 
De eek et ees tenet Mame J 
= 9, (a). AD, (er Or (a, by5 €5 x), OM (4, be, ba; C5 Y, z) 
(9) 
(tr — 5 — 5 1a, bi 5 Gay by 5 Gay ds | } 
m3) | GsX,¥,z | 
Py eee aes ms 5 


(equation continued on p. 564) 
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= AM, (cer ® (a4, by; 65 x). © (ay; a9; by, Bs" C95 Ys 2). (10) 
3. EXPANSIONS 
On making use of (5) in (8), we get 
[ES tery eee Onecian ge 
ms) | q, X,Y, 2Z 
er Ag eee eee, Cee tad 


< tnt + 
(q)r (q)s (Q)r4s 


m,n,p=0 r,s?0 


nd 


(C1)m (Co)n (€s)p (Q)m (Qn (9)p 


which (on replacing m by m +r-+ 5, nbyn-+ rand p by p + s) gives the following 
series expansions, after some simplification : 
fea OSS etl Perri + Bes Be 
om. | q; X,Y, z 
Co yp al ee eg 


-> (@)r4.5 (Dy ris (bor (D3), x7 +8 y" 24 arts g2ir?+s2)+30r 
()r (G)s (Cidr+s (C2)r (C5)s 


ee | 


r.s20 
( —:i:-; eres rts rts. fe : 
Kone aq aq’**, biq’** ; bag’ ; bq! ; 
Ll =t:—-3-—3—te,9r**, cg": cage: 
; } 
q 3 xq""*, yg’, zqr*? | mde b 


Similarly, the symbolic relations (9) and (10) yield other series expansions for basic 
triple hypergeometric functions. 


4. SPECIAL ExpANsIons 


) 
Operating upon both sides of (7) by ee (A) and then by Men (k), we get, 


(a), k 2: (b)3c,h; (d): (0), 43 (f)5(g); 


om | 
( AOC), KA): ks) 3 (e9 


9, X,Y, Zz 


_——— | 


oa) 


= 5 SS (Fras (Q~")r (Q-?), KT*8 grimentp)simtp) 
(9), (9); (kK). , 


(equation continued on Pp. 565) 


mu,p=0 rs20 
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— [(4) m4n+p [(O))mn-+n [(c))n+e (2) ]o4m [(e)]m (Cf) In [(g)]o (Am (Anse X™ Y" 2? 
(2) mtn+? [(B’)men [CC n+e [4 Norm [Cm (CF Wn (8 Vo (A) sense (Fm (Gn (Qo * 
...(12) 


Now replacing m by m + r +s, nbyn+ rand pbyp-+s in the right hand- 
side of (12), we get, after some Simplification, the following series expansion for basic 
hypergeometric funtions : 


[Oy rice (0) 7 1c), 2 310). ce), As (7); (2); 1 
] q;X%,¥,z | 
faa jena) 5 (6), ks (0) 26), SC): (s'); 


au > [(a))orzos (Cd) Iorss (+s (AD]-+2s (s+r 
Fs ((@' )}urs2s [Cb Vers s)) [Ce )rs s (a )r+2s [(C’)]s4r 


ois) 


r,s 


(h)r+5 Ch Ree \(f) [(z)], xrts y’ 2 [rts grre4 52)+3sr 
(A)orsos (K)r+s (CSD) [(g)], (Q)r (9)s 


(3) " ( ) bade (b) q , (c) q’ a gr ’ ( ) q 
x | fn d 23. 
L ( ) q ’ hq ‘ es (b ) " flighd . (c’) qrvs ( 1’) ' 


(age hg ts (fa (eas ee 
q; xq’**, yg’, zg"? |. 


megs (fos ela: 


x 


...(13) 


5. SPECIAL CASES 
Taking h = k in (13), we obtain the following peculiar type of expansion which 


has a free parameter k on the right only : 
( (a):: (0); (0); @: (23 (/)3 (8); } 
oO | | 3%, Vz | 
t (a): (6)5')3 @') 2 (2) (F95 8) j 
(k)rss [(a)]ors Qs [(c)]--s [(d)]-+2s [(e)]ras 
. S (K)or+:s [(@’)]ars2s [(c')]r+s [(4')]-+25 [(e r+ 
r,s70O 


2(r2+52)+3rs 
ena ee 





(fl (Cale (bors X77* "2? k"**g 
> [7 Jb (ds (b D.res (Gr (Q)s 
(a) get? : 3 (b) geet; (c) grt? 3 (A) git, ka: 
| 
\ 


Ahem ’ + ) qrt2s - 
(a’) qz't?s, kg? 1 (b') quits (c )qrt; (d')q 


- (0) atts. katt: (f)9" (8) 9° 3 } 
. (e)q »Kq (f)4 q; xe, yq', zq’*? 2 
(egrets (f')9' 3 (8) 4"; ) 


(14) 
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First a theorem revealing an interesting relationship between images and 
originals of a composite function of the form k (x) F [g (x)] and the function 
F (x) in two general integral transforms is established. Three other new 
theorems involving finite and infinite integral transforms are also obtained. 
By choosing the integral transforms 7; and Tz occurring in these theorems 
to be suitable integral transforms and on giving specific values to the 
functions g (x), k (x) therein, six specific theorems are obtained. Finally, an 
integral with the help of one of such theorems is evaluated. 


1, INTRODUCTION 


The present paper gives explicit expressions for the iteration of integral trans- 
forms with symmetric kernels. 


Throughout this work we use the notation 
oe) 
h(p) = Tf (x); B} = JS ke. x) f@) 7 (1.1) 
0 


for the infinite integral transform, and the notation 


hp: (b, O) = T{f ()[U & — 6) — U (x — eh; P} 


i 


i k (p, x) f (x) dx tia) 
b 


for finite integral transform. In (1.1) and (1.2) the class of functions and the domain 
of p are so prescribed that these integrals always exist. Throughout this paper the 
kernel k (p, x) is assumed to be symmetric ce. k(p, x) = k(x, p). UO - a) denotes 


the Heaviside (Unit Step) function. 
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2. MAatn THEOREMS 
Theorem 1—If 


hy [p; (6, c)] = Ty {k (x) Ay [g (x) ][U (x — 6) — U(x — cc); p} (2.1) 


and 

h, (p*) = T; {f (x); P} Be)! 
then 

hy [p; (b, c)] = 3 p $ {p, u; (B'2, 2!2)} f (uw) du (2.3) 
where 


¢ {p, u; (B''?, y"/*)} = T, {k, [ph (y*)] k [A (v2)] ye? 
x A’ (ye) (U(y — Be) — U(y — Ie); u} 
...(2.4) 


provided that the integrals involved in eqns. (2. 1)—(2.4) are absolutely convergent, p is 
independent of u, o is a non-zero real number. U (x — a) denotes the Heaviside (Unit 
Step) function, k, (p, x) and k, (p, x) are kernels of the transforms 7, and 7, respecti- 
vely. The functions k, g, f and the inverse function h = g ” are single valued, analytic, 
real on (0, o9) and g is strictly monotonic on the subinterval (6, c) such that g(b)=8 
and g (c) = y. 


Proor : If the substitution x = A (y*) is made in (2.1), then 
ce 
A, [p; (b, c)] = i k, [ph (y*)] k [A (9°)] hy [g fA (y7)}] 


x [U th (ye) — b} — U {h(y*) — oh] o ye- fy’ (y*) dy 
which easily yields after little simplification 


lle 


hy [p; (6, c)] = 6 " kK, [ph (y*)] k [A (y*)] Ay (92). po fy’ (ye) dy). ...(2.5) 
pg} g : 


Now, noting the property that each kernel is symmetric, we have k, (y,u) = k, (u, y) 
and so (2.2) yields the result 


hy (y*) = 1 ky (u, y) f (u) du. ...(2.6) 


On substituting the value of A, (y*) from (2.6) in (2.5), we get 
\/e 


Y co 
hy [p; (6, c)] = o J Ky [ph (y*)) kA) {| ke (wu, Y)f (u) du} 
B v 0 


x ye"! h’ (ys) dy, 
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Now, on interchanging the order of integrations in the above integral which is easily 
seen to be permissible under the conditions stated with the theorem, we get 


1 
a ye 


h, [p; (6, c)] = J { J k, (u, ») k, [ph (y*)] k [A (y2)] yt A’ (ye) dy} 
B g 
x f(u) du me a 
which is (2.3). 


If in Theorem 1, 6b = 0,c + co and the function g is so chosen that g(0)=8 
and g (oo) = y, it reduces to the following interesting form : 


Theorem 2—If 


hy (p) = T, {k (x) hy [g (*)); DP} ...(2.8) 
and 

h, (p*) = T2 { f (x); p} (2.9) 
then 

hi (p) = 9 F 8 (ps u; (B,D) fw) du (2.10) 
where 


¢ {p, u; (BY/%, y'!*)} = Ts {ky [ph (y9)] k [fh 0) v2 97) 
BL LY ans BUM) mm Cy 70°) | at} (2.11) 
provided that conditions of validity directly obtainable from Theorem | are satisfied. 


On the other hand, if in Theorem 1, 8 = 0,7 > co andthe function g isso 
chosen that g (b) = 0 and g (c) = ©, it takes the following form : 


Theorem 3—lf 
h, (p;(b, ©) = T, {k (x) Ay [g (IU (e — 6) — U(x — os PY 2.12) 


and 

h, (p*) = Tz {f (x); P} EAE), 
then 

hs (p:(, N= 2 5 6 (p,u) f (w) du (2.14) 
where 

¢ (p, u) = Ty {ki [ph (y2)] & [h (y?)] yr’ (y*); u} oN em BG 


provided that conditions of validity directly obtainable from Theorem | are satisfied. 
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Also, if in Theorem 1, we take b = B = 0 and let c > oe, y + oo and the 
function g is so chosen that g (0) = 0 and g (ce) = 2 [or g (0) = co and g (oo) = Oj, 
we arrive after a little simplification at the following interesting form of the theorem 
under appropriate conditions of validity easily obtainable from Theorem 1. 


Theorem 4—If 


hy (p) = T, {k (x) hy [g (x)]; Pp} (2.16) 
and 
h, (p?) = L, {f (x); P} #OA7 
then 
a) ~ 
h(p)=o f d(p,u)f(u) du ...(2.18) 
where 
o (p, u) = Ty, {k, [ph (y*)] k [h (y*)] ve"? A’ (v9); u}. .. (2.19) 
If in the above theorem we take g (x) = x, we at once arrive at a theorem given 
by Gupta*. Further if we take oc = | and g (x) = xin the above theorem, we geta 


theorem obtained by Agrawal!. 


3. Speciric THEOREMS 


We next turn to find six specific theorems from the theorems obtained earlier. 
For doing so we give specific values to 7}, T:, g(x), k (x) and co occurring therein. 
Much of the computational work involved in obtaining these specific theorems is 
straightforward, but often lenghty, so we omit the same. 


In Theorem 2, if we take o = 1, g (x) = e*, k(x) = e* and both transforms 
I’, and 7, to be the Laplace transforms, we easily get from (2.11) and the known 
result [Erdelyi*, p. 137 (4)] 


] ; 
d {p, u; (1, co)} = — ypla-bjap (2 = athe Reu > 0 
a a a 
resulting in the following specific theorem : 
Theorem 2 (a)—If 
h, (p) = L {e* hy (e#*); p} (3.1) 
and 
hy (p) = L {f (x); p} -«.(3.2) 
then 


i>) 
b 
hy (p) = + | urte-oe (2 _ z. u) f(u) “or 3.3) 
, 
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provided that the integrals involved in (3.1) to (3.3) are absolutely convergent and 
Reu > 0. 


In Theorem 3, if we take o = 1, g (x) = (x? — a?)!!2, k (x) = (x? — a*)-)4 and 
the transforms 7, and 7, to be the Laplace transform and the Hankel transform of 


order zero respectively, then from (2.15) and the known result [Erdelyi®, p. 9 (24)], we 
easily get 


d (p, u) = uw)? (u® + p?)-)? exp [—a (uv? + p*)'/"], Re p > 0, Rea > 0 
resulting in the following specific theorem : 
Theorem 3 (a)—If 
hy [p; (a, co)] = L {(x*? — a*)*!* hy [(x? — a®)'?] U(x — a); pp} ...(3.4) 
and 
hy (p) = Hot f(x); Pi ERY, 
then 


h, [p; (a, °°)] = “ ul? (u2 + p®)/? exp[— a(u? + p*)'?] f (u) du ...(3.6) 
0 
provided that the integrals involved in (3.4) to (3.6) are absolutely convergent and 


Re p> 0, Rea > 0. 


In Theorem 4, if we take « = 1, g (x) = x’, k (x) = e** and both the trans- 
forms 7; and 7, to be Laplace transforms, then from (2.19) and the known result 
[Erdelyi?, p. 147 (33)], we easily get 


$ (p, u) = $m? y-H? e suoray? Brfc [4 (p + 4) u-1/?], 
Reu > 0, | arg (p + a)| < 2/2 
resulting in the following specific theorem : 
Theorem 4 (a)—If 
h, (p) = L te7% hy (x); BP} ...(3.7) 


and ye 
ha (p) = L {f (x); P} ...(3.8) 


then 
1 2 
h, (p) = yet? f ut? hat Bete [4 (p + a) uw} f(u) du —-..3.9) 
: : 0 


provided that the integrals involved in (3.7) to (3.9) are absolutely convergent and 


Reu > 0,| arg(p + 2) | < n/2. 
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Again, in Theorem 4, if we take o = 1, g (x) = (x* + 2bx)'?, k(x) =x+b 
and the transforms 7, and T, to be the Laplace Transform and Fourier sine transform 
respectively, then from (2.19) and the known result [Erdelyi’, p. 75 (35)], we easily get 


d (p, u) = e?® pb? u (u? + p*)-) Ke [b (u? + p*)'/*], Rep > 0, Reb>0 
resulting in the following specific theorem : 


Theorem 4 (b)—If 


hy (p) = L {(x + b) hy [(x? + 2bx)!*); p} ...(3.10) 
and 

he (p) = Fy { f (x); P} iil) 
then 

hy (p) = pb? er | u(u? + p*)"? Ky [b (u? + p?)!*] f (u) du -»-(3.12) 


provided that the integrals involved in (3.10) to (3.12) are absolutely convergent and 
Rep > 0,Reb> 0. 


Also, in Theorem 4, if we takec = 1, g(x) = a(e* — 1),k (x) = (eo — 1) 
and the transforms 7, and T; to be the Laplace transform and the Stieltjes transform 
respectively, then from (2.19) and the known result [Erdelyi®, p. 217 (9)], we easily get 

,-! 
ge Petr (2 c+ 1) u 


?(Pp,u) = b 








P 
5 +1 


r(§+2) a 
A (Fthet1 2 +21- % 


|arga| <7,0 < Re(c + Do Re ( 
resulting in the following specific theorem : ; 


Theorem 4 (c)—If 


hi (p) = L {(e* — 1) hy [a (eo — 1)}: p} .-.(3.13) 
and 
2 hy (p) = S{ (x); p} (3.14) 
‘eee P (c + yr (? ars 1) f ae 
r (¢ + 2) - 


(equation continued on p. 573) 
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Pp 
ere (5 + let oe +2;1— “) £ (u) du rathoa bod 


provided that the integrals involved in (3.13) to (3.15) are absolutely convergent and 
| arga| < 7,0 < Re(c + 1) < Re(p/b + 2). 


Also, in Theorem 4, if we take o = } g(x) = (x? + 2bx)!*, k (x) = x"! and 
the transforms T, and 7, to be the Laplace transform and the Fourier cosine transform 
respectively, then from (2.19) and the known result [Erdelyi’, p. 17 (29)], we easily get 

y 1/2 
$(p,u) = 20( 5) [p+ (pt + wep (pt + uty enuestrennt, 
Rep > 0, Reb>0O, 
resulting in the following specific theorem : 
Theorem 4 (d)—If 
h, (p) = L {x7}? hy [(x? + 2bx)"!*]; p} jaet ae lO) 
and 


h, (p'!”) = F. { f (x); Pp} 3.17) 
then 


fe \le of 
mo = (FZ) 0 | tot ere te 


x exp {— b(p? + w*)'"} f (u) du my EP} 


provided that the integrals involved in (3.16) to (3.18) are absolutely convergent and 
Rep >0,Reb> 09. 


A number of other specific theorems can also be obtained but we do not record 
them here for the lack of space. 


4. APPLICATION 


On making suitable choice of f (x) in the specific theorems given above, we can 
evaluate a number of integrals. To illustrate, we evaluate an _ infinite integral by 
making use of the specific Theorem 2 (a). Thus, if we take f (x) = x** in Theorem 
2 (a), then from (3.2) and the known result [Erdelyi®, p. 133 (3)], we have 


h,(p) = 0 (s) p-’, Re Pp > 0,Res > 9. (4.1) 
Further from (4.1), (3.1) and the known result [Erdelyi’, p. 143 (1)], we easily get 
h, (p) = T (s) (p + 95 — b)-, Re p > — Re (as — b). (4.2) 


Substituting the above value of /, (p) in (3,3), we get an interesting integral 
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eo 
| yp—brse—a) ja J" ( eee f . u) du = aT (s)(p + as — b) ...(4.3) 
a 

where 

Re (p + as — b) > 0, Res > Oand Rep > 0. 
If we puta = 1,5 = | and 1 — p = @ in the above integral, we get 

co 
J. ut *-1 (a, u) du = T'(s) (s — a) ---(4.4) 
0 

where 


Re(s — a) > 0, Res > 0. 
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The classical brachistochrone problem for which the solution turns out to 
be a cycloid, is based on the assumption that during the descent the gravity 
remains uniform. So the result will be valid only for a narrow region of 
space in which the variations in gravity can be neglected. The present paper 
gives a numerical study of brachistochrones taking into account the 
variation of gravity. The shapes of these curves have also been plotted for 
Earth’s gravitational field and the times of arrivals at various points com- 
puted. The times of vertical descent and descent along straight lines have 
also been computed. Similar analysis will be useful in the field of an electric 
charge. 


1. INTRODUCTION 


The famous brachistochrone problem has a very interesting history.’ This 
problem was raised by Galileo in 1638 and its ‘solution’ though not mathematically 
correct, was a source of inspiration for much work in the field. This work was 
initiated by John Bernoulli who challanged the mathematical world in June 1696 to 
solve the following problem : 


‘Given points A and 8 ina vertical plane to find a path down which a movable 
point, by virtue of its weight proceed from A to B in the shortest possible time’. 


The above statement of the problem is followed by a paragraph in which he 
reassures his readers that the problem is very useful in mechanics and that it is not 
the straight line AB and that the solution curve is very well known to the geometers. 
He says he will show that this is so at the end of the year if no oneelse does. He 
extended this time limit to Easter of 1697. Finally, May 1697 issue of Acta Erudi- 
torum appeared with John Bernoulli's solution. This issue also contained the solution 
of his elder brother James and a brief note by Leibniz saying that he had also solved 
the problem but since his solution was similar to that of Bernoulli brothers he would 
not reproduce it. Here it may be pointed out that some accounts are available which 


indicate that Newton had also given a solution of the problem. 


* This work was submitted as an M. Phil dissertation to the Department of Mathematics, 


University of Roorkee. 
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Leibniz had suggested that the problem be called ‘Tachystoptotam’ the problem 
of swiftest descent. Bernoulli, however, prefered to name name it the ‘Brachistochrone’ 
problem (from the Greek brachistos means shortest and chronos means time). Here It 
will not be possible to describe all the methods due to two Bernoulli brothers, Leibniz 
and Newton. One may see Goldstine’ for details. The elegant method due to John 
Bernoulli was based on Fermat’s principle of least time. Today, we have mathemati- 
cal tools like calculus of variations available to us and the brachistochrone problem 
is merely a class work exercise, for which the solution is very well known to be a 


cycloid”’®, 


All the above discussion is valid only when the gravity is assumed to be uniform 
in the region in which the particle descends. The problem becomes much more 
difficult if the variation of gravity is taken into account. We have not come across 
many references dealing with this case. Smith* considers the problem through the 
Earth. He takes two points on the surface of Earth and assumes the movement of a 
particle through a tunnel dug between the two points through Earth’sinterior. It is 
well known that inside Earth the gravitational attraction is directly proportional to 
the distance from Earth’s centre. Under these conditions, as shown by Smith, the 
complete integration of Euler-Lagrange equation is possible, though through a lengthy 
and tricky sequence of substitutions. The final result is remarkably interesting. The 
brachistochrone is found to be a hypocycloid. Smith also calculates the minimum 
transit time from San Diego to San Francisco. This is found to be about 12 minutes. 
Here it may be of interest to note that the length of the tunnel along brachistochrone 
would be about 685 miles as compared with 500 miles of the straight tunnel between 
them. The time of straight tunnel is found to be about 42 minutes. 


As far as the brachistochrone problem outside Earth is concerned, we have not 
seen any papers. Here the gravity follows the inverse sqnare law and the exact 
integration of Euler-Lagrange equation is not possible. We have found a numerical 
solution to this problem in this paper. Plots of the brachistochrones have been 
obtained using CALCOMP plotter and the times of arrivals at various points com- 


puted. Comparision has been made with times taken along Straight paths and times 
taken in a vertical fall. 


2. FORMULATION OF THE PROBLEM 


Let a be the radius of Earth and let a particle move from a point P, down a 
smooth curve to another point Py as shown in Fig. 1. We have to find the brachisto- 
chrone joining P, and Py. Let O be the centre of Earth. Without any loss of 
generality, we can take the brachistochrone in the plane OP,P,. Taking OP, = r 
and OP; = ry, we should have a <'f Sr, otherwise, Keeping in view the energy 


Principle, the particle Starting from rest at P, will never reach Py. Also let the angle 
P,OP; be «. 
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Pp 


s 





Fic. 1. 


Using polar coordinates with initial line along OP,, let the position of the 
moving point at time ¢t be P(r, 6) and let its velocity there be v. Ifthe particle starts 
from rest at P, then the principle of conservation of energy gives 





ewe (=—-— 1 
mv = GMm = 3 Hi A) 


where mm is the mass of the particle, M the mass of Earth andG the gravitational 
constant. Now if g be the acceleration due to gravity at the surface of Earth i.e. 
g = MG/a’, we can write (1) as 


I ] 
ig 2 atg ( ae ...(2) 
Now using 
v = ds/dt = [r?+ (dr/d@)*]'!* d6/dr, 


where 1 is time and s the arc length from P, to P and integrating from P, to Py, the 
time taken is given by 


to = (2a*g)*” f [eee | a. 


Let us introduce the following nondimensional quantities 


...(3) 


R =rja, R, =r,/a, Rs = rela, T = t (2 g/a)'!?. ..(4) 
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Then (3) gives 


= | (oe 


R® —R,) wa BY 


Since the integrand does not contain @ explicitly one integration of the Euler-Lagrange 
equation associated with (5) leads to the following result after some simplification 


dR ( cR Ret RR, i 
dQ 


ok R R,—R 210) 


In the initial stages of motion, R will decrease with @ and so we take the negative 
sign. Then there comes a stage when dR/dg=0 and we have the point of minimum 
R. After that R starts increasing symmetrically. The point where dR/de= 0 will be 
called the critical point in the subsequent discussion. In Fig. 1, it has been denoted 
by Pc. The polar coordinates of Pc. will be denoted by (Re , 6c). Note that if Py lies 
on the part P, Pc, the study of critical point will not be important for us. But if Py 


is at P. or beyond that its position will play an important role. The critical distance 
R. can be found by solving the equation 


J (R) = cR, R° + R — R, = 0. ENT ip 


But this requires a prior knowledge of c, which cannot be determined unless (6) is 
integrated further and the two constants of integration found by satisfying the 
two boundary conditions, But some important information regarding the range of 
variation of ccan be found even without doing so. From (6), it is easy to see 
thate > 0. Also f (0) = — R, and f(R,) = cR‘,. Further df/dR = 3c R,R? ++1> 
0 in (0, R,). So there is precisely one real root of (7) in (O, R,). It is again easy 


to see that it has no other real root. This real root, i.e. Re should also satisfy the 
additional condition 


Re >, (8) 


for those brachistochrones which are not obstructéd by the Earth’s surface. This 
demands 


f (Ql) = cR, + 1-R, <0 


0 <c¢ <(R,—1)/ R, = co, say. ...(9) 


If the Earth is replaced by a point mass, it will not create any obstruction and c may 


lie anywhere between 0 and co. 


For a given value of R,, the critical value of ¢ i.e. ¢ 


o can be found from (9) and 
can be computed from (7) for various values of c. 


For example, for R, = 2.0, we 


Re. 


PROBLEM IN NONUNIFORM GRAVITY 579 


have co = 0.5 and the critical distances for some selected values of c are given in 
Table I. 


TABLE I 


Critical distances for Rs = 2.0 





, e R, c Re Cc R. c R, 
.005 1.928 .07 1.514 0.50 1.000 6.0 0.500 
01 1.869 ll 1.398 0.75 0.901 15.0 0.378 
.02 1.776 16 1.299 | Te 0.747 30.0 0.305 
.03 1.703 .23 1.202 2.0 0.689 100.0 0.208 
.05 1.595 34 1.099 3.0 0.614 1000.0 0.098 


a 


Note that for c > co = 0.5, the full brachistochrones will be available only if 
the Earth is replaced by a point mass. 


3. COMPUTATION OF BRICHISTOCHRONES 


In this section we shall compute and plot the brachistochrones through P, for 
several values of c. There will be a definite value of c for which it passes through the 
other point Py. The determination of this value ofc will be the subject matter of 
the next section. 


Since the brachistochrones can be symmetrically extended beyond the crirical 
point, it is enough to discuss the paths from the starting point to the critical points 
for various values of c. Integrating (6) from P: to P, we get 


R 
. 1 R, ale 1/2 

6= | 4 (ae ERIE) He ...(10) 
R 


Since P can be any point on the arc P,Pe, the lower limit can be varied from R- to 
R, and the corresponding value of @ can be found by evaluating the above integral 
numerically. Note that by virtue of (7), there is a singularity at the lower limit if 
it is Ree So there will be difficulty in applying the usual quadrature methods. There 
are several special methods to handle such a situation®. We have found the critical 
angle @ by splitting the integral into two parts, i.¢. 


5 in Fi ey 
Ic = i Bilas 


The second integral is proper and can be 


where € is a suitable small positive number. 
ated by 


evaluated by any numerical method. The first one can be approxim 


2 athe Re) | PR eos 4os 400) ] 
R. 1+ 3 ck, Ke’ 3 5 (11) 
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where 


1 I 3A Re | ...(12) 
ore Groseenie ort 2° 


I l 27 
Rivas Bir We Re ee 


Re 8 (R, — R.)* ae : 


3 AR. 
o 


a.{13) 
1 2R(R-R) * 4(R,—R).’ 


cR, ec ae ... (14) 
I + 3cR,Re D-2 oRe 


These formulae were used to compute points on the brachistochrones for various 
values of c for a given R,. The choice of «and the strip size for evaluating the 
second integral depends on R,. Table II gives the polar coordinates of the points on 
the brachistochrone for R, = 2.0 andc = 0.5. The results are reported at equispaced 
radial distances upto the crirical point which in this case just happens to be situated 
on the Earth’s surface. The various brachistochrones for this value of R, and several 
values of c are plotted in Fig. 2. Table IJ also gives the times of arrival at the 
various points. This will be discussed in section 5. 


TABLE II 
Coordinates and times of arrival 


(Ry = 2.0, c = 0.5) 


R 9 Time* R 6 Time* 
2.00 .0000 .0000 1.46 .0786 2.8582 
1.94 0018 9763 1.40 0993 3.0081 
1.88 0054 1.3756 1.34 1243 3.1522 
1.82 0105 1.6788 1.28 1549 3.2927 
1.76 .0170 1.9319 ie 1929 3.4325 
1.70 0251 2.1528 1.16 2417 3.5755 
1.64 0351 2.3511 1.10 ° 3077 3.7297 
1.58 0470 2.5323 1.04 4102 3.9190 
1.52 .0614 2.7004 1.00 6060 


: : 4.2100 
* 1 nondimensiona] unit of time (from (4) ) 


= (a/2 g)'/° secs, where a = 6.37x 10° m, g = 9.8 m/sec?, 


= 570 secs or 9.5 minutes. 


Note that in Fig. 2, the br 
> Co= 0.5. So they will exist 
replaced by a point mass in whic 
curves will terminate at R=2 wh 


achistochrones start intersecting Earth’s surface for 
Only upto the point of intersection unless Earth is 
h case, c may be given any value in (0, co). All the 
ere the velocity will become zero. Thus the particle 
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Fic. 2. 


will keep on moving to and fro between the two end points situated on R=2 with 
constant time period. This period will depend onc. For R,=2.0, c=0.5, the time 
period will be 16.84 nondimensional units as can be seen from Table II. There is 
one more observation to make. For small values of c, the curves are quite small in 
size and appear to be very close to a cyloid. It should indeed beso because in the 
space occupied by the curve the gravity will not vary appreciably and the problem 
will be reduced to the classical problem. However, as c is increased the curves go on 
departing more and more from the cycloidal shape as can be seen for c=1000. The 
critical points in Fig. 2 are shown circled. If necessary one can fit a suitable curve 


through these points from the numerical data. 
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4. ESTIMATION OF c WHEN Py 1s GIVEN 


In the previous section we have described the method of determining all the 
brachistochrones through a given starting point P,(R,,0). There will be one 
definite value of c for which one of these will pass through the other given point 
Py; (Ry, «). In this section, we discuss an iterative method to obtain c starting with a 
suitable initial approximation which may be obtained from Fig. 2. 


From eqn. (10) we have 


R 
iat Riek 1/2 
Of ms | (GRRE RE) dR = ® (c), say = 0. el) 
RK, 
Jf 


We have to solve this equation for c. We have applied Newton-Raphson method 
which gives 


2 V(c 
Ck41 = Ck =e, k = 0, l, hi oP: ...(16) 
where 
zo 
ee | R* (R, — R)'/2? dR 
«~ ) (RR ROR) =D 
R 
f 


A computer program was prepared to carry out the above computations for given 
R,, Re and «. For R, =2., Re = lia = 0.606, it took 30 iterations with starting 
value co = 0.3 to converge to the value 0.5, 


5. COMPUTATION OF TIME UPTO A POINT 


The time taken by a particle from P, to a point (RX, @) on the brachistochrone 
will be given by 


8 


a R? + Rk” 1/2 
“= | (a R= ) dé ...(18) 


Oo 





Where R’ = dR/dg is given by (6). Substituting this value and using the identity 


cR,R® + R-R, = cR, (R—R.) (R? + RR + K), K = (cR,)-}, UL) 
we finally get 


R 
T= VR, eee Oe 
| [(R, — R) (R= RD(R PR, R+RK)] 7° ...(20) 
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but computationally, we found the method of section 3 quite efficient numerically. 
We split the integral as follows 


R -« R 
VR,( J + f ) =I, + I,, say, if R > R.; 
R R 


> 
Ss 


and 


R +<« eal R 
e 


PRICES VET + i ae Ae a go 
R. Rit es 


Proceeding as before, if can be shown that 








2R°? el? 
I, = —_—_—— ee 
[(R, — Re) (Ri + R, Re + K) PP 
Ji +4et+ Fetow |, (21) 
Real Be 
Is sit 2 ae . 
Pike RY CYR +E) | 
[i+ 4+ Seto], (22) 
where 
2 1 B. ...(23) 
me oR By 2 
3 fi: q ‘i 
Bi Rot oa aRad) a P eg it 
a ...(24) 
+ (FR, 
p = (Re +2 R,) (R2 + R, Re + K), (25) 
q = l mt R? + R, Re 25 K), ...(26) 
...(27) 
. 2 BA. — R.) ae 3A R,{2, 
3 -2 9 2 R? ] TA 
° “3 — — Re) = ¢ ? 
B =R. + g [(R, 


c 


+( R;'- Ae) (R, — Rd". ...(28) 


584 BANI SINGH AND RAJIVE KUMAR 


The other integrals are proper and have been evaluated by using Simpson’s rule. The 
numerical values of the times in arriving at various points for Rs = 2.0andc = 0.5 
are reported in Table If in nondimensional units upto the critical point. Beyond 
that total times can be found by using symmetry of motion. It will be interesting to 
compare these times with those along straight paths and along vertical falls upto the 
same heights. This is done in the next section. 


6. TIME OF DESCENT ALONG A LINE 


Let the points P, and Py; be joined by afstraight line. It is easy to see that for 
any point (R, @) on this we have 


R, Ry sin « 


Ry sin (a—6) + R, sing * .».(29) 


R=- 
The time of descent from P, to Py along this line can be found by substituting R and 
dR/d@ obtained from this expression in (5). The resulting integral was not integrable 
in closed form. Since there is singularity at the starting point it is split into two 
parts as described above i.e. (i) from 0 toe and then e to x. The second isa proper 
integral and the Simpsons rule is applied to evaluate it. The first one is found to be 


er I 
2 Ses Ped (re )etl+ ( Sp Sot 2 ad + +O (€*) , ...(30) 


where 
R, — Ry cos @ 
Ry sin « ; (31) 


ae 


A program was prepared to compute times for various values of R,, Ry and a. 
Convergence was ensured by decreasing «. With R, = 2.0, Ry = 1, « = 0.606 
« = 0.005, the time was found to be 4.4162 units. The corresponding time along 
the brachistochrone as can be seen from Table II is 4.2100 units. As expected the 
time along the straight line is greater than along the brachistochrone. 


It may be of interest to compute time taken by a particle in a free vertical fall 


from P, to the surface of Earth. Integration of the equation of motion leads to the 
result 


° ea! ; 
T= R$" [cos (Rz"*) + R7Y*(1 —  R5*)-P], 
For R, = 2, this gives 


v2 
T= i (7 + ./6) = 3.9535 units = 37.56 minutes. 


7. CONCLUSION 


_ Although, the problem discussed here deals with brachistochrones in Earth’s 
gravitational field, the method can be easily extended to other situations where the 
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field is variable. For example, one can take the gravitational field of two or more 


bodies and find the paths of shortest time. Since it will not be possible to construct 
brachistochrones of the dimension of Earth and verify the results of this paper 
experimentally, we suggest here an experiment which can be performed in the 
laboratory. Instead of a gravitating mass we can take an electrical charge and study 
the motion of another charged particle in its electric field. Since electric forces are 
much stronger and follow the same inverse square law the brachistochrones will be of 


similar shape. Similar experiments can be performed in the field of more than one 
charge. 
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We consider the diffraction of Love waves by two parallel perfectly weak half 
planes in a layer overlying a half space. The problem is formulated in terms 
of the Wiener-Hopf equations in the transformed plane. The transmitted 
waves are then calculated using the Wiener-Hopf procedure and inverse 
transforms. 


1. INTRODUCTION 


The diffraction of seismic waves by structural dicontinuities is of considerable 
im portance in seismology because of the existence of such discontinuities in the Earth’s 
crust. Exact analytical solutions of these problems are difficult to obtain even for 
simple geometries. de Hoop? presented a method based upon the Wiener-Hopf 
technique for the solution of bedy waves by a single perfectly rigid or perfectly 
weak half plane. Kazi‘ considered the diffraction of Love-waves by perfectly rigid and 
perfectly weak half planes lying in a surface layer overlying a half space. Recently, 
Asghar and Zaman’ have considered the diffraction of Love waves by taking the 
rigid barrier to be of finite extension. 


In this paper, we set upand solve the problem of diffraction of Love waves 
normally incident on two parallel perfectly weak (crack) half planes lying in a surface 
layer and parallel to the interface between the layer and the half space. The problem is 
formulated in terms of the two Wiener-Hopf equations and can be solved by the tech- 
nique introduced by Jones®. The weak screens Separate the layer into three loosely 
coupled layers. The transmitted waves in these three regions have been calculated 
analytically. As expected on physical grounds, it has been shown that the transmitted 
wave in each region satisfies the dispersion relation of the Love waves travelling 
in a layer of uniform thickness under similar boundary condition. 





* Permanent address : Department of Mathematics, Quaid-i-Azam University, Islamabad, Pakistan. 
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2. FORMULATION OF THE PROBLEM AND THE WieNeR-Hopr Equations 


We consider the diffraction of Love waves by two parallel weak half planes (cracks) 
lying in a layer of uniform thickness H over an elastic half space. The half space 
has a rigidity «, and shear wave velocity 8, and the layer has rigidity jz. and the shear 
wave velocity 8,. The coordinate system is chosen in such. a way that the interface 
between the half space and the layered medium conicides with the xy plane, the z axis 
is directed into the half space and the two semi-infinite planes occupy z = -h,,x < 0 
and z = —h,,x <0. The free surface isz = —H. The geometry of the problem is 
shown in Fig. 1. 


Free surface Z=H 
Cree surface zcH 
V4 
ee eee eee, h2 
(inc) 
V2 
—<--— * 
Vv 
3 z=h} 
Vy vine) BP X- QXIS 
Z axis 
Fic. 1. 


Suppressing the time dependence e*”’, the incident Love waves of the Nth mode 
have the displacements : 


Vi" = A cos (o,y H) exp {—oiy Z — iK,w (x — x0)}, Z > 0 


Vi" = A cos {((Z + H) oon} exp {—iKiw (x — x0)},0 > Z > —H ..(1) 
where 
2 1/2 tke 2 yo ee 2) 
iN = (Kin ay k, ) an Son = (k; om Kiw) ! , 8, ; k, — 8, Pua | 


and Kyw is the Nth root of the Love wave dispersion equation 


5 Sol 5 1 ae oe ...(3) 


2 —— 2 S . Noto al 
tan {k, — K*H} =v (ki — Kaye” a 


corresponding to the layer thickness H. Moreover, K\y = Cas k = w/c, where Cyn 


is the phase velocity of the Love waves of the Nth mode. 


Let the total displacement field due to the presence of perfectly weak screens be 
written as 


prove! _ yiine) + V,, (j =1, 2, 3, 4) .+(4) 
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where V, represents the diffracted displacements in the jth region as shown in Fig. 1. 
The geometry of the problem leads to the following boundary conditions : 


(a) At 2 =0,00>x > - ae 
y=), 

eV, _ W, 

az BZ 


(b) For x < 0. 


oe) 
At rs = h, + 0, az : 
oy at r 
At'Z— 0, oe | 
r 37 = Asan sin (sey 51) exp 
Z= —h, +0, | 
a | — iK,w (x —- x0) ...(5b) 
A = — = ad. 
tow A ho 0, aZ | 
= H—h:i= 1,2. 
(c) Atz = — H, —co < x < co 
ov, 
a7 = 0 dal OC) 
(d) Atz = —h,x >0 
pigi pan aln exes 
(e) Atz= —h, x > 0 
V3 = V. Va _ WV 
4> aZ ie OZ" ... (Se) 


The displacements V; (j = 1; 2, 3, 4) satisfy the differential equations 


eV, , &V, rr Nao 
axe Tt age + Ks V,;=0(j = 1;/ = 2, 3, 4) ---(6,7) 





@ 
where k; = Bi’ |k, |< |k,|. The differential equations (6) and (7) together with 
the boundary conditions (Sa) through (Se) constitute the boundary value problem. 


ae With a little effort the above mixed boundary value problem can be written, in 
¢ Fourier transformed plane, in terms of the following two Wiener-Hopf equations: 


] 6h I(a, ho) 


3? sinho.h I (a, h,) se (a, rae h;) + Ke anaes (a, my _— = 4 ; 
2 2 


(equation continued on p. 589) 
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{Viele —h) + Vi_ (a, —h) | = h (Ve. (e, —h) —Vp- («, —h,)} 
...(8) 


» 
— 


} 
ooh ie (a, —Ai) + Vy_ + (a, hy) } 








62 sinh o,h 
— 1 oA __ sinh o 5, 3, 
o2 sinho,h 36,5, 38, 
< 2 (a, — h.) + Vy_(4, —h,) = {V3_ (a,—hy) —Vq- (a,hg)} 
..(9) 
where 
a? = (a? —k*i)'/?; I (a, hy) = 92 sin 9 hy + vo; cosh oh 


h=h, —h,, 5; H— hi. 
script ‘+’ are analytic in the domain 


In equations (8) and (9) the functions with sub 
Im (« + k,) > 0 and those with subscript ‘—’ are analytic in the domain Im 


(a > sae k,) < 0. 
3. DETERMINATION OF THE WIBNER-HopF SOLUTION 
the 


the Wiener-Hopf equations (8) and (9) can be obtained by 
Omitting the details of calculation, we can write 


The solution of 
usual procedure outlined by Noble®. 
the diffracted field in various regions as : 

3, cosh ¢,Z — vo, sinh o,Z H, (a) (< oh i) 
a—k 


V¥3(s, 2) = — o, sinh 6, h, + voi cosh 6, hy Ts («) 
iA oN | T. (Ky) sin G26 91 
| Ns ©) — Fam Ge Kin) PON 2) Rant ha) He (Kio 
_..(10) 


gin oy 3s [Ps (2) — Po (Kin) | 


(3 ak aie 





_ cosh ¢ (Z + H) H, (a) 
Vs(a,Z) = sinh os H T, (4) \a — kp 
iA oon exp (ikin XO) ¥ T.. (Kyw) sin (¢2n) 9 
~ (a — J (K\w + k,) H+ (Kin) 


|». (a) — = 
— gin (on 52) [P+ (x) — Ps (K,y)] \] 
( a + ke ) * 


a—k, 


cosh o, (Z+h,) Hs (*) 
sinh o,f Y, (4) 


(equation continued on Pp. 590) 
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. i 8 
iA N exp (iK,n x0) Le (K,n) sin GoN Op 
| 0, (a}= = : a SEa (Kisastia) A, (K,n) 
— sin (o,~ 5,) [W, («) — Wy (K,»)] }]. 7G) 





Ve lt 2a = 


cosh 5,(Z + H) Hy, (a) (: as ae 
a-—-—k 


sinh 5, 3, Y, (a) = i 


iA Gay exp (iKiw x0) _ [¥¢ (Kiw) sin (ay 81) 
[ 0, 0 ar ea Raa ee 


— sin (%2n 3) [W, (a) — Wy, (K,w)] i] Patio) 
In obtaining these results, we have used the factorization 
Snh Cet =H, (2). H_ (a) 
ooh 
— * = T, (x). T- (2) ...(13) 
soa Ge) a ty = Ys @).¥(@) 


which have been explicitly obtained in the Appendix (A). Also, the splitting technique 
of Noble® is used to write the following additive decompositions (the explicity forms 
are given in Appendix B). 


Fay (x, —h,) 
‘ae ke) Alea «) + N_(a), 
(a + k,) A, (a) T.(a) ~ N+ (#) + N_ (a) 


(@ + ki) H, (2) T- (a) ~ 7+) TP 6 (14 
(% + ks) Hy (a) T- (a) ~ P+ (2) + P- (a) (14) 


Vix (%, Aj) : 
(«+ ky) A, @) ¥o(a) = 0+ @) + O- (a) 


] 
(« + kz) Ay (a) Y- (<) = W, («) + W_ (a). 


4. THE TRANSMITTED WAVES 

We determine the transmitt 

half planes in the layer. This ca 
for each region separately, 


ed waves in the three regions that are formed by the 
n be done by taking inverse transforms. We do this 


(a) The region — fhyheo Boer Nste < 0 
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The Fourier inversion formula gives 


€+00 
] 
v, (x, Z) = a | exp (—i/ax) V, («, z) dx veal 45) 


c- oO 


where V,, (a, z) is given (10) and Im (k,) > c >— Im(k,). For x <0; -—h, <2 <9, 
we can close the contour in the upper half plane. Lapwood® showed that the contri- 
butions to the surface waves come through the poles. The branch point contributions 
give rise to body waves which are of no interest to us for the present study. The 
integrand in eqn. (15) has simple poles at « = K,y and at the zeros of / («, h,) located 
in the upper half plane. Using the relations (2) and (3), the contribution ¥,,; from the 
pole at « = K,n can be written as 


¥os, (X, Z) = — A exp {—iK,w (x — x0} cos {o2n (Z + H)} ...(16) 


which cancels exactly the incident Love wave as is to be expected. Let Kam 
(m = 1, 2, 3...) denote the zeros of J («, hy) in the upper half plane. Then using 


Cim = (xe — ki re, Som = (k, x Ko 


we can write 


io.2) 


feo Bete aH) Kye ; (Ce -1). Serax 
i cae >» Poa a rar gh (Som AD \ 7 xp (—#Az, 


m=1 


K,,, + ko 1/2 Fa. (Kom) [ (K. ) = iA Gon exp (iKin x0) 
| os k, Ty (Kom) pies 20 (Kim — Kyy) 


: is (K,n) sin oon Oy a ne $. (P. (Kom) —Px (K,w)] }]. 
x li. (K,w) (Ky + ks) sin o,v 5, [P, (K: ) + (ALN 
Le 
In (17), we have used 
| oF. h, 
dI (a, h,) | I ae =e (18) 


Com _ 
da E=K Kom cos (Sam hy) (3 rat, l ) 


gm 


‘ se At 
where Com = is the phase velocity of the Love type waves ol the mth mode 


a i i = ich i uivalent to 
may be noted that (17) requires the relations 7 (K,) = 9 which is eq 


Sim (19) 


=s Vv . 
tan 7, Ay Se. 


structure 
Since (19) is the relation for the propagation Love type waves es pat reer 
consisting of a semi-infinite solid of rigidity 4 covered by a surface 1a) 
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thickness /1, and rigidity 1, v2,. is therefore a Love propagating in the geometry as 
show in Fig. 1. 

(b) The region —-h<z< —h3z<0 

The transmitted wave v, (x, z) in this region is determined by applying inversion 
formula to eqn. (11). The contour of integration is closed in the upper half plane. 


The integrand has simple poles at « = K,y and at the zeros of sinh oc, éA lying in the 


upper half plane. The contribution v3,,; (x, z) arising from the pole at « = Ky, 


cancels the incident wave »,"*in this region. The contribution from the poles at 


nn? 1/2 a 
a = ip,, where P, = [ Fy tas Kk ] , denoted by v3,, is given by 





7H 
< Det bP py fe ZtD 
Vay. = +/2ni SS exp (Pax) 2 > | HM. (iPn) S wie 


n=0 


iA Gov exp (iKin x0) ia (K,y) sin (sew 5,) 
2n iPn _— Kin (Ky “+ ky») H, (Kin) 








— sin don 5, [(P, (ipn) — P,. (K,w)] } 





nr 
salad Saas aI _ iA. egw exp (iK,y x0) 
Ys (iPn) * 2m (inn — Ky) 


ns (K,w) sin (oon 5) : q 7 . * 
x eee, — sin (6,0 3,) (W, (iP,) —W, (Kw) TT 
...(20) 


Note that in deriving eqn. (20), he have used the relation sinh 3, h = 0, which is the 
dispersion relation for waves in an infinite strip of uniform thickness § and rigidity p, 
with weak upper and lower surfaces at z = —h, and z — — h, respectively. 

(c) Theregion -H<z< —h,» x <0 

To find the transmitted waves in this range, we apply inversion formula to eqn. 
(12). Closing the contour of integration in the upper half plane, we find that the in- 
tegrand has simple poles at « = K;y and at zeros of sinh o,8, that lie in the upper half 


plane. The contribution, »,,,, from « = Kin exactly cancels the incident wave ven in this 


n? = 





region. The poles arising from the zeros of sinh o. 8. are « — ip, = { 3 — 








52 
n= 1,2, 3,.... These poles give rise to the contribution 
oo ao ° 
ip, + kz 12 14 (ip!) 
Vaso = 4/2ni s exp (p‘ x) cos — (Z + H) [ 
n=0 ; ip’ — k, Y, (ip, ) 


(equation continued on p. 593) 
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Bre ro fe cap Ux) {¥. Kin) sin oa by 
/ 2n (ip’, — Kin) (K\w+k.) H, (Ky) 


— sin ozv 3, (W. (ip) Wy (K,w)) }]. (21) 


The dispersion equation satisfied by 4,2 corresponds to the relation sinh oc, 5, = 0 


chat; ‘ n® a? 112 ee 
atisa = | kz — —3z | » which is the dispersion relation for waves in an infinite 


Strip of uniform thickness 5, and rigidity », with free upper surface and free/weak 
lower surface respectively. 
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nN = 


Pom 


APPENDIX A 


The factorization of the function, involved in eqn. (8) and (9), has been fully 
described by Sato’. We only quote the results here. 


(a) Let us write 


co 
a |e i} = Aw (AN) 
2 
nei 
where 
4 2 a i. 3, 
5, 3 = [1 — 3 in ~i 2-2, B=. ww(A.2) 


Also H (a) = H,. H_ (a), where 
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Hs (a) = Tl pf, 8, F in By + exp (F iad + X (a)}. (A.3) 
n=) 
If 
foe: » 5 BS ~ p 
% (a) = — ia o| t-e- log ~ +>} then Hz (a) ~} a]? 


as | a | co 
in appropriate half planes. Hence 


sinh o, 6 
o, 5 





= H, (x). H- (2). ...(A.4) 


(b) If + Kimand + Kem (m = 1, 2,...) denote the zeros of / (a, A.) and I (a, hx) 
respectively, we can write 


A 
”» 
a? — K- 


ren) It (52) 











where 


G,(«) = 1 (a, h)| Me — K2) 


G, (a) = I (a, h,)/ It (a? _ K?,) 


and G;,, («) has no zeros. 


Let 
O@= GS = 0,@).0-@) 
then 
log (= > | tog sO) da, OES) 
: 
Tak} = I. (a). O(a) = Ty (a) T_ (a) —...(A.5) 
where : 


Ts (x) = Ox (a) <j 


. (A, 
Lae (A.6) 
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(c) The parallel calculation from (Al) to (A4) for ah nS St af and sah e¢% 





lead to > Shad ay ds 

sinh a, 4, Gy 5, 

3,8, sinha, 8, ~ ‘+ (*)- ¥- (a) iA) 

where 

sinh o, ¢ 

Sake TIO) % (0) 

sinh c, 5 

aa = Yee (a) - Ya- (a) 
and 

_/ ¥, + (a) 
Fae as.(3): 


APPENDIX B 


Following the general decomposition theorem (Noble®), the explicit respresenta- 
tions Ns («), Ps («), O+1 (a) and Wx (a) are given by 


id+ 0 
ic+ o V! , —Ah,)d 
Wrrey crc w 5 Sense tok : ...(B.1) 
+ 2ni (E + k,) H, (8) T- (8) (& — +) 
ic—oo 
hi dé 
| Sei pepe ee ..(B.2 
Pelham = i \ere H, (&)T-(&) (& — 4) ae 
vi, (&, —A) dé 
a I ee ye fs Fe 
Os() = + 577 lernHoroe= 9 ee 
dé ...(B.4) 


1 ens ee aoe 
Wa (e) = mal + k,) H.(@) YE — » 


where —Im (k,) < ¢ < Im(z)<d<Im (ki). 


The integrals in (B1), (B2), (B3) and (B4) can be calculated by the contour inte- 
gration method. Let us consider 


; 1c4.00 Vi, (8) 48 
v.0=p7 | €FHHOC-OTO 
oo 


fc- 


...(B.5) 


—y, cos 6, where |g| <7. Thus, eqn. (BS) 


= 


For « lying in(—%, n) we can write « = 
can be written as 
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[ove) A 
[e+c0 Viv (61 7 h,) HG aa Kon) dé 





l 
N+ (—% cos @) = Tai | “(E+k,) H, (E) je (£) (E+n cosé) ...(B.6) 


{c- 


4 
r 
s 





where T_ (&) = Q_ (&) Il [ 


m=l1 


of T_ (€) in (B6), we obtain 


= | and Q_ (&) has no zeros. Putting the value 


hand K, 
E ao K;,, 
{c+00 ge (E; a hy) nG TS Kam) dé 


] 


co A 
fe~00 (E-+k») H4(&) O_(&) n (E — Ki,,) (+7 cos 6)" 
...(B.7) 
Closing the line of integration by a semi-circle in the upper half plane the poles 


captured are & = — 7 cos 0, & = Kim (m = 1, 2...). Thus 


co A A 
EZ Ve. (Ku —hy) MH (Ku — Kem) 
N,(—7 cos @) = _— 


A 


ma (Ky + k,) Hy (Ky) O4 (Ky) it (Ki = BS) 








Vi, (—1 cos 6, —h,) 
(Kz —7 cos 6) A, (—yc0s 6) T- (—7 cosé) 
_.(B.8) 


te 


The other integrals can be evaluated similarly. 
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Let S be a surface and S be the spherical image or inverse surface or parallel 
surface of S. By knowing the relations of the first and second fundamental 
magnitudes of S and S, we study the properties of streamlines of Sin S. 


1. PRELIMINARIES 


For steady flow the system of equations are given by® 
(oy (Ch)'= 
(ii) yx Vx H=V X Ga (Vv X H)) 


Git) pV DV me DH-¥( ptm) +07 


(iv) PV -vy=V-V-T)+Y -«(VN+(V x 
x [pe ve (VY X H) — pe V X H| 


(v) H=«V (1.1) 
where 

p = density, T = stress tensor, 

H = magnetic field vector, T = Temperature, 

V = velocity, je = magnetic permeability, 

vy = magnetic viscosity, h = stagnation enthalpy, 

p = pressure . = heat flux, 


a = scalar function. 
Using (1.1) (iii) in (1.1) (vy) and ultimately using this in (1.1) (iv) taking account 


of the vector identity 


a’ 
(a . 9) a= curla x a + grad (4 ) 
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and (1.1) (v) we get 
y2 
P(V - yh) =VvV .- |e (curl V x V + grad (=-)) — pe (curl (@V) xX oV 


2 2 2 pe 
+ grad (“So ))+ v(p +H 2s ) } tye 


+V-(VT) + (9 xX A) (evn (y x H)— peV x #). 








veh iee) 
Using (1.1) (v) we have 
Vx H=0 Estiaa) 
curl (ec V) x (2 VV) = «curl V x V ...(1.4) 
For two dimensional flow 
BT ety = Tyr C11 + Tz Crp + Tor €or + Too Co * a tye 
Using (1.3), (1.4) and (1.5) in (1.2) and properties of vector analysis we get as 


: ' V2 2 2 
P(V - vh)=phb grad (-) — ne ¥ + grad (SY ) 





2 2 


V 
+V (p+ me “5 Jtv-@vD 


+ He va (YX H)? + Tye, + Tys C19 + Trg 22 + To) ey1. 


re saa 420) 


V = V(E, 4) e. mt? 


In natural i.e. streamline co-ordinates with 81 (5 ) dE and gp (E, n) dy as the 
components of a vector element of are length we have the following system of 
equations (1.8) 


var 1 @g 
(i) div (PV) = 36 3 (g, pV) 

Laces a Sy 01, ae a oll 2a 

Zecw kr 06 cea) gare )e 
wa a? V2 12 22 3 
(iii) grad ( 5 )- (5 Jet 3 ra] (e% Je 


(iv) grad h = 1 OA 1 oh 
2, 


(ii) grad 

















OB See nee 
ih ails 1 @ 272 
(W) grad (p+ pe * )= = oe ( a’ V Lig? 
2 Pam 3 i Baan) 935 ohm 
a? V2 
x (P+ ne 2 ) es 
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ENP ee ee ie 
g, Fahd Th» gz, ay ey 
7 ‘ oT ; ne 
« grad T = — = « oF 
& 0& a 22 on % 


2(- oT  \, O[« a 
ae\g, & * )+ mem s:) | 


(viii) curl H = y xX H=vy X (2V) = grad« x V+ x curl V 


Gil) div @ ered T)—"— 
81 8s 





1 ACe 1 @ 
gradx« = —-—>e Sa Bee 

pie ey On. 
grad Ie ea Ae 

22 On 
| 7) 
curl V = — ante V 
2, 2 on (g, V) es 


V ou ¢. ra) 
yx H=-— A ee oes 
a ae ee 


1 @ 


a= — Fg, on xV)es 
, | a 2 
1X x HYP = —>wZa wis A 
(ix) (V ) oi gi ( a @ Pa v)) (1.8) 
Using (1.7) and the properties of stress tensor in (1.5) we have 
2 eV wee. Ww ag 
e = errheao >, = — —— — = ae — 20 
a wo ee a3 22. On ( 21 ) anc ey 21 82 0& 


and 


2 

2 2vn eV \ Pat ae WeLga th 
myare 2(ore BE) +m HCR(H) 
> “LS Z dias ae) oe + 2vH 92 \ an\ gi 


81 
{-1 
( 2vy Vi 8g, \ 8s 
Sat g (Se $ ee Fe) wash ee) 
£182 e 81 82 0 of 


a unit vector in the direction to the velo- 
dicular to the velocity but in the plane of 
1 to the plane of flow such that ¢1, &» 
f this (2, 7) net is of the form 


Here V is the magnitude of velocity, ¢; 
city, €, a unit vector in the direction perpen 
flow, és a unit vector in the direction of norma 
and e; from a right handed system. The metric 0 


ds? = gi (E, ) 4% + 8 (E, ») dn” 
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where g, and g, satisfy the Gauss equation 


ue ba $: ) 6 (2 08 ) 
—» (| — = | + | — -} = 0. 
“te & 6m \8. ay 


Using (1.8) (ii) to (x) in (1.6) we have 


PV eh _ 1 2(é oT )+ Ei ges )] 
s 8 ~ ges Laels, & )* Or les on * 


ana (= ( xv )) (> 
gig? \ay \% Near: ae a) 


nr aed a? V2 YA 2vH OV OV 
we? 3 | 2 Narersa tae: . ) 3 


+ 

















8 8 o& 
g? ( a Ge VAY 2 
ea ACC) Far 
"et Ven\ ag, J) + 2g 8s 
2 wH 08> 3 O82 
a «Dies» 
( p+ aS (1.10) 
From (1.1) (i) and (1.8) (i) we have 
6 
aE (g,P V) = wCeehs) 
and this implies 
82 PV = f(a) 
implies 
Pea db he 
82P 


Using (1.11) in (1.10) and cancelling the term I/g, g> we have 


AO — ele aee)+ Ree) 
tone Cin P ae) + oar we ee) 
GO (EO)+ Bett (Le) 
+ AE os Bf te) ap 
+ (rs te 2 (£9) 2 (29). 


Ba 
- 








8 
ea 
SNH Si a 








GRP) 
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We prescribe the streamlines to be straight lines. We assume that they are not 
parallel but envelope a curve C. We now take the tangents to the curve and their or- 
thogonal tracjectories the involute of c as a sytem of curvilinear co-ordinates. The 
aban element of arc length ds in this or-thogonal curvilinear co-ordinate system is 
given by 


ds? = dt + (E — o(n))' dy? (1.13) 


where cs = oc (y) denotes the arc length of the curve and 4 the angle subtended by the 
tangent line with x-axis and § is a parameter constant along the involute of this co- 
ordinate system, the co-ordinate curves § = constant are the orthogonal trajectories 
(streamlines) are the involutes of the curve C and 4 = constant (the streamlines) are 
its tangent lines. 


Hence by’ 


gi(6&,n) = 1, g. = §& — ¢(y). ...(1.14) 


2. Tue Stupy of STREAMLINES OF MGD ELOw OF A SuRFACE S IN THE 
SPHERICAL IMAGE S 


The study of spherical image is a spherical representation of a point, differential 
cross-section or partial configuration of the surface on a unit sphere whose centre may 
be considered as origin. The movement of a point on the curve on a given surface has 
a corresponding relational movement on the sphere. It can be mentioned in this con- 
nection that one of the recent developments in nuclear magnetic resonance spectros- 
copy is to consider the spectrum by a surface in 3-dimensional space in constrast to 
the higher to hold notion that it is a graph in terms of absorption intensity as a 
function of frequency. This is one case where in the concept of spherical representa- 
tion can be inculcated to study the spherically nuclear magnetic resonance spectral 


properties. 
Let S be the spherical image of a surface S. Let n be the unit normal at the 


point P on the surface S, the point Q whose position vector is n is said to correspond 
to Por to be the image of P. Clearly Q lies on the unit sphere S. The position vector 


R of 2 on S is given by 


R=n. at2i) 
Differentiating (2.1) w.r.t. 6 and 4 we have 
R, <n, = W7? ((FM — GL) + (FL — EM) r.] sar(202) 
R, =n, = W? ((FN —- GM) r, + (FM — EN)r), 8: = - ; 
éR (2.3) 


R,= an, 
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where r is the position vector of P and E, F,G and L, M, N are the first and second 
fundamental magnitudes on the surface S and (&, %) is a curvilinear co-ordinate system 
on the surface S. 


Let e, /, g denote the first fundamental magnitudes for the spherical image 5: 
They are defined by 


e= R, : R,, f= R, Rs and g = Rs . Rg. ...(2.4) 
Using (2.2), (2.3) in (2.4) we have 


e = W-* (EM? — 2FLM + GL*) (25) 
f = W-? (EMN — FM? — FLM + GLM) .«.(2.6) 
g = W-? (EN? — 2FMN + GM?), X27) 


Let the streamlines in the surface S be straight lines. Using (1.14) in (2.5), (2.6) 
and (2.7) we have 


e = (& — a (r))-? [M2 + (E — o())? LJ ae a.) 
f= (& — o (0)? [MN + (E — o (y))? LM] ...(2.9) 
g = (6 — o(n))?[N? + (E — a (n))? M?]. .-.(2.10) 


The image curves of § = constant and » = constant of S will be orthogonal in 
S provided f = 0 i.e. M = 0. Hence the image curves of § = constant and 7= cons- 
tant will be orthogonal in S if & and y are Jines of curvature because F = 0 and M = 0 
is the necessary and suficient condition for the curves — = constant and 7 = constant 
to be lines of curvature. In this case if k, and ko are the principal curvatures, 
these are given by 


L=k,E, N=kG ax(2.01) 
where k, and k, are functions of &, 7. Hence using (1.14) we have 
L=k,M=0,N=k3(& —o(n)). - aan) 
Using (2.12) in (2.8), (2.9) and (2.10) we have 


e= ki, f=0,g =k? (E — a(n). »- (2.13) 


Remark 2.1: If the lines of curvatures § = constant and 7 = constant are 
parametric curve then 


ek, 


ak, kk 
on 


of Rah Oke ] k Sarr k, . 
E E, and aie 2 Cie aes G, sal ecad 


ty | — 


where EF, and G, 


denote the partial differentiation of E and G w.r.t. 1 and & respec- 
tively. We h 


ave seen that the flow in the surface S is planar. 
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Similarly we can show that the image of the flow of S in the spherical image § 
is planar by using (2.13) and (2.14) and 


I a 1 avg g 1 éaéVve 
K= — ——— — —_—-| -———. = 0 
/ eg ( =r (- Je oe ) a ot (= on )) 
Take e = gi (&, 1) = ki and g = gi (&, n) = k? (& — o (n))? and those imply 
81 (8 0) = ky, 82 (8 0) = ke (E —  (y))?. of: [5) 


Theorem 2.1—If the streamlines in MGD in the surface S are straight lines but 
not parallel and envelope a curve C and the tangent lines to C i.e. 7 = constant and 
the involute of the curve C i.e. § = constant are taken as orthogonal curvilinear co- 
ordinate then the image streamlines of S in the sphere S are concentric circles. 





Proor : Using (2.15) in (1.12) we have 


oh a Kk 
p(y = (4 e-em 2) 


é k hs Me YH 
‘ ee ean on he + Ee eH eG) 


x A (eS elee eo) 


f(r) a? f (1) AU es 
oe Sore =( p2 k2 (E = ()* )+ ky (F (3 sind (n)) 


d f () ‘ 2 f (1) 
les ol erase oo) + pk, (§ — ¢ (m)) 


x| = p+ aes | « ( ts 2) + ky) 


Pers 2a) 


¢lsne=a ) 


Simplifying the above equation we get 




















eT « k, 
é é « k, OT hs ks amet “Paes 5) 2 
fo) =e - %( eo a ao 

‘ a2 T 

O42 k Oe a 

+ ala )E=o & te 


(equation continued on p. 604) 
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x oO (n) k, eT He YH e + ( k, a f (n) ) 1 


+E —c}' on | Kika —e)LOn\ Phe 


k, «' f (0) * f 3% (n) 6 
Te Re = ae ee eaGiics), 


(ear > PRES |- 4% 


[Ge (ae gee eee 
+ EE ot) os (pie eae eee 


2f (n) ‘ 2va f (a) 
FoR eo) | Mt 


F—o 


S 




















+ IG — | —p + 3H atten 
~ Pk, Ena} Ol seGe je 


] 
Pk2(& — a)? iS 
Multiply the above equation by (& — o(n))® we get 





ae — - ‘al K k, er k 
ES o)? elaltage( yt | eee 


hy a eT a : 
+ ——* (- — o/s BP eee? ee RA, k, oT 
BSP ea a Sra ft ky ye 





ck (Eo) @T . .o' (yk aT 





re ee ge (HELO) & - 0) + ee | 


jo nae 


rar) © - o) — aaa E — 0} | 
_ Hef * (n) 


op Lae ( peaE)  - 9) ~ pee 0) ] 


(equation continued on p. 605) 
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de 8 (f(a) f(n) 3’ (n) 7? 
aa > itr 8 ra ‘Ei 


2f(y)yF _ 2vu f (0) [( ok. 
et pee ae Ct) 





ok. 
x | ae (G-= ) +-.k, (E — °)] + 2k, (£ — 6)? f (2) 


x | - PAG ia eaten a 


x (ER - a) x le -~ (se) | 
... (2.16) 


Equation (2.16) to hold identically we must have & = o(n), Usingthis in (2.16) 
we have 


neva KPa (eV S20) » 2a ky fPmo*(). _ 9, 


k, ke p* k2 k, p? k? k2- ees eae 


Simplifying (2.17) we have 


" ; 2 
tle (n) 0°? (n) vu (He ky he + —) ==) 0) 
1 
implies 
f? (n) (4) vx = 9. 
Since f (1) = 0 and vy # O we must have o’ (n) = 9. 


This shows that streamlines are concentric circles. 


3. Tue Stupy OF STREAMLINES OF MGD FLow OF A SURFACE S IN THE 
INVERSE SURFACE S 


As the name inverse surface itself is suggestive of is the surface corresponding to the 
given surface which has the same shape, size when viewed from a crorresponding point 
on the given surface. In other words concavity and convexity characteristic of a given 
surface assume the converse characteristics on an inverse surface. An analysis of inverse 
surface with that of the given surface helps us to study a surface at different orienta- 
tions Gaseous flow, incompressible, compressible all encounter such sort of flow spots 
simultaneously and consecutively. It is true that a perfectly geometrically congruent 
isometric natured conformal charactered matched surface and inverse surface is 
difficult to come across in practice. An approximate evaluation of the pattern of flow 
can be made if such an assumption were made that a crest and crevice, a mountain 
and precipice are taken to be set of given surface and inverse surface. With the 
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establishment of the new fundamental magnitudes for the inverse surface and know- 
ledge of intrinsic relationship between the fundamental magnitudes which determine 
surface and its inverse surface, it is possible to maintain the same tempo of flow. 


Let S be the inverse surface of S. Let the center of inversion be taken as the 
origin. If c is the radius of inversion, the position vector 7 of a point on the inverse 
surface S corresponding to the point r on S has the direction of r and magnitude c?/r?. 
It is given by © 


2 
fda My, .G.1) 


Differentiating (3.1) w.rt. § and y we have 








. OF _ c 2c? 
i; = ae os = i a Ne Si Pee (53 
Sak Se as 2c" 
eat ema ThE girs ees. Pa ts} 


Using (3.2), (3.3) the first fundamental magnitudes, E, F, G of S are given by 


4=T7, > ry = pe E, R = Pr) . Fy = re F, 
ee cé -_ e 
=r °'h = 7A G, and W= ae W, cosh ace 


By assuming the streamlines in the surface S to be streamlines we establish that 
the streamlines in the inverse surface S are also Streamlines by using the relations of 
first and second fundamental magnitudes of S and S, 


Using (1.14) in (3.4) we have 





ci = c! ro o 
E = pe? F=0,6 = ek (6 —o(y))? and W? = = (6 — o())2. 
Se & ie 


, We have already seen that the fiow in S is planar the Gauss equation correspon- 
ing to the fundamental magnitudes of S is given by 


vara e+ iso) 


VO. On (3.0) 


Z 4 
Take gi (, n) = E = <andg?&,) —G = & (E — 6 (1)? 3.7) 
implies 


=e CS Cae 
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Lemma 3.1—If flow in the surfa i 
I : ce S is planar then the i i 
surface S of S is either elliptic or hyperbolic. pee aot ae Ee 


PROOF : Using (3.5) in (3.6) we have 
R~ — eae Lal al (E-° @) 
+ lees oey a ())] 
Sat aga e\ a )) 
+ oy (g=sor &(e))]~ - eee 
«Lal BS) ht alee 
ogee 


- (sew) 


ee _ 2n 2-5 (a) 2(—E—o()) - 
. GEO r a? ee ry pee as 











— 2 





ea Saeed 2rs eee ek ee 
r(& — s (m)) Fs ree 5 (1) eta | 


= 2r3 ; e — o n) 
mad ae ey |" +(G@-—s (4) ru — oy ry 
ED a. re - re ] 
— — a(n) (—E — ¢())’ r(6—o (ns 
Thus K is either positive or negative. Hence the image flow in the inverse surface S of 
S is either elliptic or hyperbolic by 


a 


Theorem 3.1—If the streamlines in MGD in the surface S are straight lines but 
not parallel and envelope a curve C and the tangent lines to C i.e. 4 = constant and 
the involute of the curve Cie. — = constant are taken as orthogonal curvilinear Co- 
ordinates then the image streamlines of S in the inverse surface S of S are either con- 
centric circles or the magnetic permeability of the streamline is constant. 


Proor : Using (3.7) in (1.12) we have 





oh é oT Oo A OT ; r‘ lke VH 
f(a) aa ae 2 oD (<e ) sa eas a Oy )+ c* (& — 9) 


(equation continued on p. 608) 
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«Le lre= a] +a" aelrere =e) 


- 








a He f (4) e( sal ddl 6) se ) . 2vy r° ae 

2p GE \ pect (& — a)? c (& —*a) rs 
he YP 4a@re’—l i 26-9) , 
Pc? (eeccey) cr FES s) r3 ; 


«| et eee Pete = ee e- On) 
« alse) ae) 


Simplifying the above equation we get 








BB ong OF ge ce ee Be eT 
f(a) ee oe (5 =o) by a 





wie Be BES. Mee RT, < er 
E—o oy a | Eo a + Ene Oy 

rh Be Vi f(n) a af'() af (n)o’ 
+ ope, : 4 FO) 0’ () 





a 








P(E —a) dy PA(S)— @) P (5 —is) 


_ 4f(y) oP Z| f3 (a) [ 4r° ry 2rs 
Evo Fr 2c8 LP2(E — 0)? ~ PFE — op 


ee es roe | - SC teen 2 rt a? 
PG ain ae 2pct po(E— a)?  pt(E— op 


6 2h, op 2 ar! 2vH 
peop  * peop e Solem eee 


cot We 4r® ry f(y) f(n)r* ep 


“P(E—o) P(& — a) P2(— — 6) dy 


6 para ré} 2¢ (n) ei Ge, 
= — a)? ree ( — s) 


i peg afi “. ~ s) r)] 

















r} 


Pe (ne) 


42 -E= 9) = arr; r? ap r2 
7100 | et PE—sc) PE —c) +ceom | 


(equation continued on Pp. 609) 
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609 
ee out maa ae aay Fe 
3 rea ) , ...(3.9), 
Multiplying (3.9) by (E — P)°, we get 
& — oF f(m) = E — 0)’ Legs (E — 0) oe OF & —o)t 7 
46-7 0G — op TF 





OT | rt we vu fn) a. 
+ (E —- a)4 Ox? + mee eae | | Gs) (> an 


af’(n) _ «f(n) 6p a f(n) 0’ (n) 7? 
fee mote 5 


f? (a) 3)? 2r* ry a ees 
+ FM @- of | AE - 0) -4 


r3 ¢ op = He f* (n) 
— piG — 9) oe Pct 


4 4 ra? @p 6a Ee 
«(Te (atr- Bae az) ) 
2M | Ef (a) + Arf (a) 


- Lo F)E-g4 “LO | 
+ FO) @ — o[1- =8 =), | 
ee pe (1-28 9) 


2 y( — o)* 2rry —o 
+ LOG =o) A 9) 





(& — s)? 





r? «2 
oo |t 


2 2r? vi f(y) 
> oP — — || —p -«)* + —*4 
as 56. % Al P(&—9) od 


ee ee 
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Equation (3.10) holds identically on the curve — =o (») and finally we must 
have 


r* pe YH a? f® (n) 9? (4) iP =i cbe Tal UU Je Pama .G.11) 
c! p* r P 


This implies that 
va fo" oP (4 + 2)=0. 


Since r 2, vy ~ 0, f(q) ~ O and p # 0 
we must have o° = 0 or pe = —2c*. 


Thus the image streamlines of S in the inverse surface S are concentric circles or 
the magnetic permeability of streamline is constant and it is equal to —2c*. 


4. Tue Stupy OF STREAMLINES OF MGD FLow oF A SURFACE S IN THE 
PARALLEL SURFACE S 


A surface S which is at a constant distance along the normal from another sur- 
face S is said to be parallel to S. As the constant distance may be chosen arbitrarily 
the number of such parallel surfaces is infinite. If r is the current point on the surface 
S,n the unit normal to that surface and c the constant distance, the corresponding 
point 7 on the parallel surface S is given by 


7=r+ten. .(4.1) 


Let the lines of curvature on S be taken as parametric curves. If a, 8 are the 
principal radii of curvature at the point r on S. The corresponding point A on the first 
sheet of the centre surface is 


rF=r+an. ...(4.2) 


Now since A is the centre of curvature for 7 = constant, it follows that 7 and « 
are constant for one differentiation w.r.t &. Thus 


r, + an, = 0. re tial 
Similarly 
ro + Bn, = 0. ..(4.4) 


Since the lines of curvature on S are taken as parametric curves, so that F = 0 
and M = 0. Then by virtue of (4.3) and (4.4) we have from (4.2) 


(c — a) 
r =< sisieas ry 


1="m ten, = = ...(4.5) 


ro = fo + ech, = — Jes ~ PB) r 


8 a: .. (4.6) 
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The magnitudes of the first order for the parallel surface are given by 


1 A Cora | ad 
\  ...(4.7) 


ro (ap renee ala | 


A : 

pet eae streamlines on the surface S to be straight lines but not parallel. By 

: wing the relations between the fundamental magnitudes of S and S we study the 
ow of MGD of surface S on the parallel surface S. 


Using (1.14) in (4. 7) we have 

















E= (<=), 0. Go= (— Fy) (E — o(n))? (4.8) 


where «, $ are principal radii and k, «= 1, k 
2 a=, ky 8 = 1 and* ka, k faci 
curvatures. 6 are principal 


Let E = gi (&. ») (cka — 1)*, F = 0, 
G = (ck, — 1)? (& —  (y))’. 
Then g, and g, for S are given by 
gi (&, n) = (cka — 1), 82 (8 0) = (cho — I) (6 —  (M). ...(4.9) 


We know that the flow in S is planar. The Gauss equation corresponding to the 
fundamental magnitudes of S, is given by (3.6). 


Using (4.9) in (3.6) we have 


aot «(4-5 a (ck — 1) (6 — 6) 
£18206 \(cka — 1) Ot ) 


CA(e=oeso* 
= raat —a(¢ Fe (6 — 20 ) + (che — ») 


+2( I eH) 
én \ (cko— 1) (6 — oly) om 


K= 


Use (1.14) in (2.14) we get K=0. 


Hence the flow in the parallel surface S of S is planar. 

e surface S are straight lines but 
nt lines to C Le. 4 = constant and 
ken as orthogonal curvilinear CO- 


Theorem 4.1—If the streamlines in MGD in th 
pe a curve C and the tange 


not parallel and envelo 
= constant are ta 


the involute of the curve Cielté 
=k), kyo=kz. 





*Note : From now onwards treat Ka 
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ordinates. Then the image streamlines of S in the parallel surface S of S are either 
concentric circles or the magnetic permeability of the streamline is constant. 


Proor : Using (4.9) in (1.12) we have 


K 


aha eT 
f(a) oe. - (sea ae hs < rape 





r) 

a rel Ge iT) © — 6 (7)) z a (oh 5 ) 
Me VH é 2 (= (cka — 1) af(n) )] 

* Ce — 1) (cke — 1)(E —c — 6 (n)) On \ P(ckn — 1) (& — o) 

f(m) 4 BSE ee ik oe? He f(n) 9 
ae: Peat — o)* (cky — or] “2B BE 
~ ( . a® f* (1) ) 2vH (cka — 1) 

P* (cky — 1)*(— — a)? (cks — 1) (& — a) 


é f() 2 2f (n) 
79 Pam eres Ci be Tt (n 
Eat — I) (cko — 1) (E - 5) + P (cks — 1) (— — o) 


ee FAG), 
x[-p Tp (ck, — N(cky — IP Ea op 7 (ck — 1) é-<) | 


x & (chs —1)(& — )) + 2(cky — 11) (E — a) 


x[- a Ree hoe elm hezs)] 
xz ® geal = nie ...(4.10) 


Using (4.9) and (2.14) in (4.10) we get 





ch _ (cky — 1)(— — a (m)) & 8T” Es 
f(y) ee NORE = ANG (y ie x(cks (—E—s) eT 
oe (che 1) SO de? ee eens 





«(cky — 1) @T | cx(&—c) & ka—ky 
(cka— 1) @ (cka— 1) @ (E—o) 














_ te (cky — 1)(& — 6) #k, aT c(cka-1)) @T 
(ck, — 1)? Uta? (cko — 1) (& — c) 67? 





__(cka — 1) o« aT cx(cka — 1) _0T 2k, 


(cke — 1)(& — 0) &q Gy ~ (che — 1)® (Eo) an an 


(equation continued on p. 613) 
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k (cka — 1) 0’ (4) or Ite VH 
* (chy — 1) — 0)! Oy) * (ha — D(cke —1) © — ©) 
. al f (4) cke — | « (cka—1) o' (4) OT 
éen\ pp J} (cko— 1)(E— 4) — (cky—1) E—0)* On 
be YH (ee cka — | 
Tk. -- 1) (cks — 1)(E — 9) p (ck»—1) (E—3) 
af(n)e 0k, _c(cka — l)af(n) Oke 

















P (cky — 1)(E— 0) Om ~~ P (chy — 1)? E— 9) Om 


4 a’ (4) (cka — 1) ee ey Qvn (cka — 1)8 
P (cka — 1) (E — )? | + (cks — 1) E — 9) 





‘ EE (ee) _ a TeaENy (— —o) 


a J(n) 3’ (a) cf (n) ok, 
(cka— 1) (cko—1) (E—0)2 P(cka—1)? (ck»—1) (6-9) on 
cf (9) __ 24m) 
~ P (cka — 1) (cke — “1 (& — 9) a = | + p(ck»—1) ‘p(cko— 1) (E—9) 
2vu f(n) 5) Oke 
e |- as (cka — 1) (ck» — 1)? (& — 9)? {e@ - ) Oe 





Pky Dt] [ee — 0) RP + (eke - 1] 


a (itn) 
+ Ack» — 6 ~ |e (AP ae DES 
of (1) bk, __c(cko — I af(n)_ eke 


p(ck» — 1)* (& — 9) 0 ~~ P(cko — 1)? (6 — 9) “on 
Qvuy (cka — 1)* 


o'(n) (cka — 1) af (4)? ioesanie Be 
+ Dick — NE — of | * che — DE - 9% 





f(a) 
x o(@) eee ER 
f(n) 9’ (n) a A ee 
+ ack=1) (che—1) G= 2) (CRI) (eho 1 B= 8) 


ok 2f (”) 2vH f(a) : 
* On | * p He =o} |-? + p(cke—1)* (6-9) 


x (cke — 1) + 2cko — VE - of a(S) 


(equation continued on Pp. 614) 
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, l ae cf (4 (ka = ky) 
(ck, —)E—«) ~ Pek, — IPE — a? 


f(a) 
~ P(ck, — IE — OEE Pe —-1) 


xf —f(n) cf (n) (ka—ko) 
P(cks — 1) (E — 5)? P(ck» — 1)? (E — a)? 


(9) | 
A a (e p HS cee = 1) @-—9) i]. ...(4.11) 
Multiplying (4.11) by (€ — o (y))°, then we have 








: # _ (cks — 1) 6 Ox oT 
(& — o(y))° f(4) Se (cka — 1) & — CO" ce GE 


« (cky _ Ib — id be oT 
ta (Ckaeet) amen es +6 — oF 





0 
_ r(cko =) gg Ake eT | «(cke — 1) 
(cka — 1®* *° GE 8k icky — 1) 
pf UR Retr, ee (okaeek) dx OT 
GU HT eat) oy an 
_ cx (cha = 1) — 0)! aF ke 
(ck, — 1)? On Oy 
r « (cka — 1) 4’ (n) (E — o)8 eT 
(cky — 1) oH 
He VH = (He cka—1 
# Vek. ae — 1) (cky — 1)| ay p a5 Co) 
_ e(cke — vate) ek, 


pk — 1, © ~ Os, 


3° (yn) (cka — 1) % (4) (cha — 1) af (y) 2vH (cka — 1)3 
13 P(ck» — 1) mer ayer Barer 7s ToS 5 by = ~~ (cke— 1) 


ae (A) eed ae 
: = ( p oe — 1) (ck» — 1) 
f(a) 6'(n) f(1)(—) aks}? 
+ Picks = l)tckp==1)” Pick ia 
P(cka — 1) (cks — 1) P(ckKa—}) (eho —1* a 
2f (1) (& — 0)? W 


(equation continued on p, 615) 
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xX (cka—1) + 2(ck —1) (E—s)? l e(4) (6 — o) 








e&\ pe )( cke—1) 
_ of (4) (cka — 1) 
P(ck» — 1) |[- p(&— «+ kee i) 
(& —e) 2 (f(y) )(cka — 1 
‘ Fes — |} 2 (4) = poe oP}. 
...(4.12) 


The above equation (4.12 i i ¢ 
ater hoc q (4.12) holds identically on the curve ¢ = o(7). Finally we 


be YH (= (n) (cka — 1) f)V" 4 2vn (cka — 1)* 


(cka — 1) (coko — 1) P(cke — 1) (cke —1) 


__ fa @) 2 
j: (ar. = 1) (cko es 5) =e 


The above equation implies 





weve ——__ (a (n)\* (cha — It f2 (0), vw (ck — I) 
(cka — 1) (cko — 1) Pp? (ck, — 1)? a Be ~ a 
f? (0) (a (a)? 
: p2 (cka — 1)* (ck» — 1) ae 
implies 
ka — 1 
Wen aal (n) 0°? (0) vit (He + 2) = 9. 
. Ka aoe l 

Since f (n) + 0, vw + 0, p #0 and a ice 
we must have 

a’? (n) (ue + 2) = 9. (4.13) 


This equation (4.13) implies 
o’ (n) = Oorpr +2 = 9. 


Hence the image streamlines of S in S are concentric circles or the magnetic 


permeability of streamline is constant and it is equal to —2. 
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